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THE ANALYSIS OF TURBULENT FLOW BY HOT WIRE SIGNALS OF POOR QUALftY 

M. Bartenwerfer 

1, Most technically important flows of gases or liquids are 
turbulent flows. As a rule, these flows are stationary in a time- 
average, so that the quantities of state of the flow medium are 
composed of time-independent averages and stochastic fluctuations, 
for instance: 

Pressure p(x,t) = p(j<)+p' ( x,t ) with ‘ p'=o, 

Velocity u(x,t) = p(x)+u'(i.t) with ;u'=o . 

Flow fields can only be described phenomenologically and 
theoretically when the speed is measurable. In contrast to laminar 
flows, for turbulent flows, not only is the average speed u to be 
determined, but the fluctuations must be characterized, e.g. by 
specific time averages. Depending on the type of theoretical 
description, these can be correlations ^^7"^ (i,j = i,2,3) as important 
quantities of the Reynolds stress tensor puTuT or even correlations 
of higher order, e.g. (i, 

Many measurement methods suitable for laminar stationary flows, 
for instance, for the Prandtl static tube, are not suitable for 
measurements in turbulent flows, since the correct, average speed 
is not measured and since average values of the fluctuations cannot 
be measured at all. 

Now it is a known fact that a warm body in a cooler, flowing 
gas cools off faster the "stronger the wind" and the greater the 
temperature difference between the body and the gas. According to 
this principle, hot-wire anemometers measure speeds and their 
fluctuations in flowing media. The potential measurement of 
temperature fluctuations will not be discussed here, rather — 
except for the small, turbulent temperature fluctuations — a constant 
temperature shall be assumed at the measurement location. 

A thin wire tied between two prongs is heated electrically. /V 

In order to keep the temperature of the wire constant over time 
during the initially laminar flow, more heat per time unit must be 
generated, the greater the rate of flow of the moving gas. The 
heat removal is affected significantly by the speed component 
perpendicular to the -wire, and only a little by the component parallel 
to the wire — if it is not neglected entirely. At any rate, the heat 
flow grows uniformly with these two components, if we neglect any 
possible influence of the prongs and/or the probe shaft. 

Through suitable electronic circuits — summaries e.g. in Brad- 
shaw [1], Strickert [2], Comte-Bellot [3] — it ispossible to keep 
the wire temperature or heating current of the wire constant in a 
chronologically changing speed, or in a turbulent flow, and a 
starting voltage of the anemometer will increase uniformly with the 
heat flow, or uniformly with the speed. Through comparison measure- 
ments, advantages and disadvantages of the constant-temperature or 
constant-current anemometer were derived, e.g. by Helland and van 
Atta [4]. At present, the anemometers with constant wire temperature 
seem to be preferred. 
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As already mentioned, as a rule the hot-wire is cooled much 
more by a perpendicular flow than by a parallel one. If the wire 
diameter d ratio to the effective wire length 1 is sufficiently 
small (d|l < 1|200), then for many probes or anemometers, a good 
approximation is: 


£ ?■ = FKr/ ( uj- + u/ ; 


( 1 ) 


with E 5=^ 

Eo“ 
FktO = 

ui 

irn 


operating voltage of the anemometer, 
operating voltage of the anemometer at u = 0, 
monotonous function of the argument, FKt(O) = 0 
velocity component perpendicular to the hot wire 
velocity component parallel to the hot wire (see fig. 


1 ) 



Velocity-independent constant k*«i, k*»o .04 for d/i»i/200. 


For other probes, using an accurate observation of the heat 
transfer over a large velocity range or when noting the prong-in- 
fluence, in the literature equations other than (1) are given, e.g. 
by Hinze [5], Webster [6], Corrsin [7], Champagne et al . [8], 

Davies and Bruun [9], Friehe and Schwarz [ 10 J . 


/lO 


Surprisingly, regarding the used equation (1), the parallel 
component Um cools the wire much less (k^ << i) than the perpendicular 
component u| but it acts with the same exponent 2. Simmons and 
Bailey .[llT found in 1927 that the heat loss of an electrically- 
heated filament increases linearly with the rate of flow in a 
parallel-flowing stream, but in a perpendicular stream, it increases 
linearly with the square-root of the speed. Accordingly, if an^ 
additive composition of these different cooling effects is possible, 
instead of (1), a relation 

- Eo*" •» Fki { + fe* w#'*' ) 

would be expected. Probably, the relation (1) still proved useful 
because in most cases, the tangential cooling represents only a corr- 
ection. In particular, in modern hot-wire probes, the tangential 
cooling can be neglected as a rule. 

If eq. (1) is true, then by using an ahalogous linearization 
step or by using a digital computer, even a proportional relation 
can be expected between the output signal E^ and the square 
of the *feffective" velocity Ugff.^=/[,Y+k^u, ;■■■:!• : 

£* » • ( u/ + k^ - u„* ) (2) 


K = velocity-independent constant, K>0 (unit: Vg/m) . 

This linearization simplifies the analysis of hot-wire signals 
in turbulent flows . 


The experimental validation of the functional relation (1) , the 
setting of the linearization and the determination of the constants 
k and K is usually performed as a calibration in a laminar or slightly 
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turbulent flow. In measurements in turbulent flows, it is assumed 
that the equations (1) or (2) also apply for instantaneous values 
of the speed and of the signals, even though the speed in a turbu- 
lent flow can change very quickly. This assumption presumes that 
the electronic control of the anemometer can follow the fluctuations 
at sufficient speed — in order to keep the wire temperature or the 
heating current constant — and also that the cooling of the wire can 
be described as quasi-stationary . 

Thus, by using one point as an example, it is clear that the 
hot-wire measuring method is not as simple in practice as indicated 
by the summary presentation of measurement methods. Actually, this 
method is based on decades of testing, use and development. A 
bibliography by Freymuth [12] contains about 1300 publications on 
the subject of thermo-anemometry and this illustrates the scope of 
research work in this area in the period from 1900 to 1978. 

Like all measuring methods, hot-wire measuring has its own 
sources of error and limitations. In the course of time, many 
problems have been corrected by refinement of probes and electronic 
equipment, but rising expectations of measuring accuracy and util- 
ity are also noted. If we presume the validity of a relation (1) 
or (2), we limit ourselves to incompressible, turbulent flows of 
constant temperature (except for the turbulent temperature fluctua- 
tions), or to low flow speeds (flow mach-number «i ) . Vagt [13] 
reported on the state of the art of hot-wire probes for just this 
case and also discussed the occurring sources of error in detail. 
Thus, a summary of possible errors will suffice here for illustra- 
ting the prerequisites for further investigations. 


2. Errors 

As for all probe measurements, it must first be assumed that 
the probe does not affect the flow. In general, this assumption 
is correct when both the sensor (hot wire) and its mount (prong, 
shaft) are sufficiently small or thin. The meaning of "sufficiently 
small" in this case, the best shape of the prongs and shaft and 
the size of the remaining error, must all be determined by experi- 
ment . 

A sufficiently small spatial size of the hot wire is also 
a prerequisite for a local measurement. The speed along the wire 
must be viewed as constant. A velocity gradient perpendicular to 
the wire may not be allowed to deform it. In multi-wire probes, the 
speed in the entire measured volume containing all wires, must be the 
same. On the other hand, interaction between the wires must be 
excluded. The optimum hot-wire lengths and spacing will depend on 
the turbulence structure of the measured flow field. For incompress- 
ible flows, the diameter of the measured volume should not exceed 
ca.*l mm (Hinze [5]). Figure 2 shows the outline of a modern X-wire 
probe which meets these requirements. 

*Translator ' s note: ca . = about. 
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Even for optimized prongs and shafts, they will still disturb 
the axial symmetry of the stretched wire. It may then happen that 
the sensitivity of the wire will depend on the angle of incidence 
(of the probe) 7 . In simpler cases, for 7 <|- (see fig. 1), we have 
the following expression instead of ( 2 ): 

£2 » *■ 1^'' ) ; ^ 2 ) 

^ ^ I ~ ^ \ '1 . 

Neglecting this directional characteristic leads to errors when h?^l. 

Due to the flow, the wire or prongs can start to vibrate; natur- 
ally this will falsify the signal compared to a fixed wire. Prong 
vibrations are usually recognized because they occur at a certain 
speed-independent Eigen frequency. Through suitable prongs, an 
appropriate prestretching of the hot-wire and through a useful 
alignment of the probe, these vibrations can be prevented, but in 
case of doubt, a control measurement is needed. 

When taking measurements in boundary layers, or when the average 
speed depends on the location, and when as a rule turbulence quan- 
tities are to be determined, the hot-wire must be positioned with 
great accuracy. Both its local position and its direction with 
respect to a specified coordinate system must be known with sufficient 
accuracy. Turbulence quantities can be erroneous by an inaccurate 
angular setting, especially when signals of a hot-wire in various 
orientations are needed for their calculation. 

If the probe is set near a wall, then wall influences can impede 
the measurement. First, the probe in connection with the wall, can 
disturb the flow so much that it will be changed significantly. 

This must be prevented by a suitable probe design and a favorable 
alignment of the shaft. Secondly, for heat-conducting walls, the 
heat is removed from the wire by differing physical laws than in a 
free field of flow. If this effect is ignored, very large errors 
can be made in the signal analysis. Wall influences cannot be 
determined theoretically today. It remains to be explained what 
minimum distance from the wall is needed to permit neglecting of 
the wall influences on the hot-wire signal (see Vagt [13] and works 
cited there). 

The characteristic line of an anemometer e=£(ui> Uti >•••) depends 
on many parameters which must remain constant during calibration and 
measurement. One of the critical quantities is the ratio of wire 
temperature to the surrounding temperature. The electronic control 

itself can be not quite stable under some circumstances — it may 
drift. As long as these changes are not completely understood, only 
a calibration before and after the measurement is of any use. 

It was already mentioned that the hot-wire anemometer is 
'statically calibrated" as a rule, and that the resulting heat-trans- 
fer law also applies for instationary flows. By comparitfve "dyna- 
mic .calibration" using soundwaves or shaking table, this assump- 
tion has been verified by experiment; see Bechert [14], Bremhorst 
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and Gilmore [15], but only for low frequencies. 

To measure turbulence quantities it is recommended to use a 

linearization step, since the analytic relation E(u) is then simple. 

A linearization can increase the noise level and the drift 
of the entire system; frequent calibration is thus indespensible . 

The linearization is as a rule only accurate for a certain speed 
range, whereby the low velocities are critical . In high-turbulence 
flows, the instantaneous values of speed can lie outside the 
"calibration range", even when the average speed lies within in. 

An ideal linearization is assumed here. 

In the analysis of hot-wire signals of turbulent flows, low 
turbulence intensity is presumed throughput and a series expansion 
of the signal E by the small relative fluctuations uj/U is used. 

It is clear that through this restrictive condition, a systematic 
error increasing with increasing turbulence can occur in the speeds 
and their fluctuations. 

This long and yet incomplete list of sources of error illustra- 
tes that only very carefully performed hot-wire measurements will 
give reasonable results. For a more accurate discussion of the 
individual errors, refer to Vagt [13] and to the works cited there. 

A part of the error is probe-related and development-induced, e.g. 
the prong influence; another part is of a fundamental nature. But /15 

over-all, the number of possible errors is very large, and many of 
them cannot be estimated realistically. These facts impede a ^ 

quantitative error computation and thus one is not available, ■ — 

Experiments have shown that several errors can be of considerable j 
magnitude and can var;y greatly. This can be the reason why 
systematic errors have not been taken into account very much. 

In this paper, the analysis of hot-wire signals is discussed 
with regard to the named systemmatic errors which follow 
from the series expansion of the speed fluctuations. Possibilities 
for an exact signal evaluation are also sought. Furthermore, we eval- 
uate the size of the error in the relative fluctuations 
resulting from the use of one- or two-wire probes after the 
unavoidable series expansion. All other measurement errors presented 
above are ignored. It will be shown that the use of one— and two- 
wire probes is subject to fundamentally lower limits due to high 
turbulence than has been previously assumed. 

It is also assumed that the frequently used equation: 

;E ^ - K* • ( • u,,^) ; (2) 

K,k=speed-independent constants , 
and that it applies for*: 

’"One could also begin from eq. (3), but essentially the same results 
are obtained under the standard conditions k*«i and |h*-i)<<i . 
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-the instantaneous values of starting voltage and flow speeds, 

-any incident-flow direction to the probe, 

-all speeds below the maximum speed of speed of sound). 

Physically this means that ideal, linearized, hot-wire probes 
must be studied which 

-are not subject to prong influences, wall influences and undergo 
no prong or wire vibrations, 

-operate at constant temperature, 

-are found in an incompressible, turbulent or laminar flow of con- 
stant temperature. 

These extensive prerequisites must be met in order to obtain 
quantitative results. The discussed, system errors may be just 
as important in more general cases. 


3 . Exact Analyses 
Let: 

£ = + ( 2 *) 

In a rectangular, right-handed oriented coordinate system, 

Kz’ is will now describe the direction of a hot wire by 
the unit vector d in polar coordinates with respect to the xj^-axis: 


5 = cosoijc^ + sina cos9 

+ sina sinO 2 I 3 • 

0 i a i -2 > 

«Tt < 9 S ?( . 



The polar angle a is the angle between the x^-axis and the wire. 

It lies in a range which will describe all possible wire direc 

tions. 


The speed u is broken down into components: 



U “ 


£1 U, X, 
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With : 


^ X 

-C u, ** 

A 



"* “ 

u • 4 


• cos^ + ■ Smoc • J'/w ^ 

T3 

<-oS 


follows for 

the 

quotient e=E/K from (2'): 

' 3 



xz ^ 

^ ^Xn, • > 

tSB 



= /f 

+ ( • eos^Kf , 



- A 

+ 1 k’- - -i) ■ sm*' !' ■ Ct>s^ ^ > 



<=■ ✓f 

+ i yt) ■ j -«'* & > 


*-f2 


{ < • Si*t t< ■ toj ^ ^ 

j 

3X3 

«r 

(. - a) • sin^<< ■ C4>J 0 X/,, 0 ^ 



- 

— a) • o< * << • JrVf ) 


(4) 


X»» 


— 

rt>i 


The other goal of an exact analysis must be to obtain ai system 
bf equations 

3 

£ 


_ a, ( ot: ) ■ u. ■ 


/ » <... 


// 


(5) 


through selection of several directions 
which can be resolved by the three quantities: 


O' 4 u ' 






If this resolution were to succeed, then we would have complete 
information about the velocity, namely the time history profile 
of orthogonal velocity components. From these, one could derive 
any time-averages, especially the compone nt s of average velocity, 
UIj U2j Uo, the turbulence quantities ujuj, u^7[, 

and also time averages of any order. 


If we select the usual orientations for X-wires, a=|; e^o.Tu |,/ 2 , 
-n/ 2 ! n/A, -3-/i/Ai 3 ti/a, -w/a and denote the signals in this sequence 
as e. to eg, then we obtain the following system of equations 
(derivation of these equations and explanation of matrix notation 
in Appendix A) . 

- k* H (U^ = (ei -<>/^ - (4 ^ ~ j 

(g.Z) + ^ 


(C.S) ('f -e^+ep/ 2 -{£ ^ 




A-fk^ 

A 

z 


'(t'. 


f X z ' 

J|C+ 

‘i 


A-fk^ 

d. 

T 

‘t 




ej- + e,^ 
k 




3 + k^ 


(Ui+u'f 

' 




e 

S ~J 



1 <5’ J 
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The system of equations (6.4) is not solvable because the equations 
are not linearly independent. In appendix B it is shown that a 
complete resolution of the linear equation system (5) by the formal 
six unknowns will not succeed as long as « is held constant 

and only 0 is varied. Equations (6.1) to (6.3) give exact formulas 

which will be needed again later. 

Now if we choose six suitable orientations for the hot-wire, 

e.g. ;(«,6) = (Ti/G,n)! (i:/3,7i): (7r/6,tr/2); (n/3, -n/?); f(n/6,ir/4 ) j (Tt/3, -3 ti/4 ) 

we can actually resolve the system of equations (5): 


(7,-fJ - 

n.3) -h up ^ ■= 

i 

bjL <’<i, • 



( * u; ) ( (4 + ^ 

4 

i 



ly.sr) 

( Ug +- j * ( O 3 + uj ) 

k 

bz3 

> 

Lf.t) 

L C -F uj) 

f V ^ 

1 

e. * . 


The coefficients are given explicitely for k=0 in 

appendix C. We see there that qll ;'b^^(a., 0 .)^o , so that in equa- 
tions (7.1) to (7.3), all six ejf actually occur. Assuming that 
at any time is true*, then by time averaging the primary /20 

sought quantities can be found: 


{t.4) - 

< f. S) 

«« ( 

£ 
i » W 


, I = -r, 

- 

IS. C) 

up ■=■ 

C 


up , r» ^,2,3^ 

IS. 7) 



1 

r 


r 

(S.») 



w 

£ 

9i)~^ - 


U.3 ) 


b* * 

£. 

1 

bf^ Ip(.\ 9i) cp — 

y^y, ■ 


If the average speed components Ui , U 23 U 3 are known, then 
equations (8.4) to (8.9) provide all elements of the Reynolds 
stress tensor. Only the time averages e.^ of the hot-wire 
signals are needed, but no correlations ^ 7 ^’ stationary 

signals, these averages do not have to be measured simultaneously, 
rather, they can be determined in sequence. But the determination 
of the average speeds is difficult since it seems as if the 
signals of six hot-wires are needed simultaneously. On the other 
hand, equations (7.4) to (7.6) are entirely dispensable if from 
(7.1) to 7.3) it follows: 

6 , ^ » /2 

SO that from this, ^'1 (i = i.2,3) and “1”^ (i.n-=.u2.3) can be derived directly. 

The system of equations (7.1) to (7.6) is thus overestimated (likewise 
the system of equations (6.1) to (6.4)). One can thus hope that an accurate 
analysis is possible with fewer than six simultaneous hot-wire 
signals. 


*The sign problem is discussed in section 3. 2. 1.1. 
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3.1 Laminar Flow 

If the velocity field is time-constant then so also 

are the hot-wire signals. Equations (6.1) to (6.3) in this case 
give : 

~ = te“ - < ) , 

t /f -k'-; C4. Uj - ( e/ + e/; 13 . - ^ ^ 

( /( -K*) Ci, '» ( 


/21 


Thus one can easily find a coordinate system by using the criteria: 


*■ 


with U 2 = = 0 and = U>0. 


From (6.4) it follows: 


^1* 



or 



(9) 


For a single hot-wire which is sloped by the angle “ to the 
speed, in general: 

^ ± ( 10 ) 

J1V7 « + J{^ co.r'^fV 


3.2 Turbulent Flow 

Initially it seemed that a hot-wire probe with six wires 
would be-needed for a precise calculation of the average speeds. 

But actually we do not have 6 unknowns, but only three: 

= U^+u^ . U2 « U2 + u^ . Uj = + 

This suggests that in principle, three different hot-wire orienta- 
tions will suffice for an exact analysis. The general case of a 
three-wire probe with three different hot-wire directions does not 
have a closed solution, but this is not an insurmountable obstacle, 
since we can find many special cases in which the desired resolution 
by velocity components is possible: namely, probes whose three hot- 
wires are perpendicular to each other or those which lie in a plane. 
These cases will be discussed below. 


3.2.1 Orthogonal Three-Wire Probes 

Now imagine three hot-wires in the direction of the orthogonal 
coordinate axes h and their signals denoted as e^, e 2 , e^. The 

system of equations (5) is simplified into the linear system: 
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1 / 

4 

" 1 




'f h 





( i/j +1.,' j 



The solution to this jsystem ot equations gives: 


\ c t/, 1.; / 


< 5 + k*)( 


K*) 


e/ + e/ 4- Cj* ^ 

tt +- 



( 11 ) 


( 12 ) 


Let us assume for the moment that for all times t: 

UjU) ^ S: O ( l ^ Z,3) 


(13) 


thus, the velocity vector i! 
u^,u 2 ,u 2 >o . Then we have 


(u^.u^.uj) does not leave the octanden 

+ Cj >-•• Cj + e/ , 


and furthermore: 





S'/ 


Iw •“ 'f, ^ t ^ > 


(15.1) 

(15.2) 

(15.3) 


This result is formally quite simple. It was derived without 
approximations 5 the condition (13) is an easier condition than 
that needed for a series expansion |uj'|«u. , but also means a 

restriction on turbulent intensity. Since the equation (14) is a 
linear relation between the fluctuations in speed and those in the 
hot-wire signals, eq. (15.2) does not apply precisely to frequency 
components . 

The turbulence intensity up and the correlation *J'u; U^m) 
were computed with respect to a coordinate system which is defined 
by the three orthogonal hot-wires. This coordinate system may not 
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coincide with the flow-line-related coordinate system (y^^, y 2 j y^- 
axis = tangent, normal, binormal of the flow line). Now 
if a component U;]^ disappears, then the corresponding condition (13) 
is violated even for low turbulence. It is most useful to place 
the main flow direction into the symmetry axis <ii+i 2 +f 3 ) of the 
coordinate intersection, so that _ is ^ ^ 


and thus j^i = ^2 = ^3 

Now if the speed fluctuations are to be determined with respect to 
a flow-line or a fixed coordinate System, they have to be transformed 
again after measurement (see section 8 ). 


3. 2. 1.1. The Sign Problem, Probe Symmetry 


Now when is u^(t)>n at all times? What are the consequences 
if an existing change in sign is ignored? How can one find the 
correct sign? 


Viewed statistically, there can be sign changes in every tur- 
bulent flow in the jvelocity components. When the average velocity com- 
ponents > 0 , then only in case of low turbulence is the time 
interval with Uj_ < 0 negligible. In general, from equation (12) 
it follows: 




and not: 




If changes sign in the time history profile, then this is /2k 

not noticed by a single hot-wire or by an orthogonal three-wire 
probe; both "see” a speed fluctuation where negative velocity values 
are handled as positive values of the same amount: 



Key: 1-fluctuation "seen" by the hot wire 
2 -actual velocity fluctuation 
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Using the example of a simple periodic step function, one can 
see how the average value and the signal scattering are invalidated 
by this reflection: 



— 

u (t) 



!u (t)i 




; 



— 

— 

— 




1 

■*- scheinbor 

1 


J 



U 2 

2 

wirklich 





'^wirklich 






t — ► 


Key: 1-apparent 2-actual 


In general it can be proven that the absolute value of the signal 
|u(t)| has a higher average value and a smaller fluctuation than 
the actual signal u(t), when it assumes positive and negative values. 


The system of equations: 

3 

e.-^ =“ ( Ot.-, 0^i - u - 

X^nt tx /i 


(5) 


is invariant against sign reversal of . Thus, there 

can be no unique solution to the system of equations by u, rather 
with u a -u is also a solution of eq. (5). When the qua"3rattc 
forms ” -'>1 rn diagonal, i.e. for l?^m, then (5) is 

also invariant to sign reversal for each individual component. 

With ji-Uu, then also are solutions of (5), so 

that there are 2^=8 solutions. This is the case for the orthogonal 
three-wire probe. The ambiguity of the solution thus has to do 
with the symmetry level of the probe: A very symmetric probe gives 

less information than a less-symmetric probe. For this reason, 
it is recommended to examine non-orthogonal , asymmetric probes 
(sec. 3.2.2), although the solution of (5) is more difficult for 
l“.hese probes than for orthogonal probes. Naturally, the ambiguity 
of solutions is limited since the speed components are continuous 
functions of time. 


Also note that the ambiguity of the functional relation between 
hot-wire signals and the speed components plays a role when using 
one- and two-wire probes when they are used in flows of high tur- 
bulence intensity. In theoretical investigations by Tutu and 
Chevray [16] and by Bradbury [ 17 ] , whose results are presented in 
section 6, this so-called rectification error was estimated and 
recognized as important in high turbulence. 
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3. 2. 1.2 Alignment of Orthogonal Three-Wire Probes 

Theoretically, only the wire directions are specified, but 
not their precise local position. The three hot wires must be 
close together in order to resolve even small turbulence structures; 
on the other hand, they must be spaced far enough apart so that they 
do not interfere with each other, e.g. one wire should not lie in 
the wake of another — related to the average flow. Figures 3 and 4 
show possible probe designs. These probes are invariant against 
rotations by 120° or 240® about the i^^-axis (probe axis). The 
wire directions are perpendicular to each other, the angle a between 
the wire directions and the probe axis is 54.7° (precisely: 
la z. arc CDS (3)"*^^= arc sin ((2/3) ‘''^)the Same angle as between the ledges 
of a cube and the spatial diagonal. As long as the angle a between 
the instantaneous flow direction and thh probe axis is sufficiently 
small, at no time will one wire lie in the wake of another: 

^ ^ ( ( j'i + «.i (^. 1 '}^ ] 


for the 'trian gular" probe (fig. 3), e.g. A s 35.5° for l'=0.5; 

A < a = 5 ^. 7 ° for the "tripod probe" (fig. 4), regardless of 1 ' . 


The linear transformation between the base vectors x. (1=^ 
1,2,3) of the orthogonal coordinate system specified by the hot 

wires and the base vectors £i (1 = 1 . 2 , 3) with as probe axis, is 
defined by: 




j 'If-tJ 

0 i/vr \ 


> 1 ^ 


'tf-fF 

jo 

¥z 

Ui 


i ^ //r 


1 y» / 


for the "triangular probe" or by 


> A 1 


1 


0 ' 


( V. 1 

A 

- 

^/■fT 


-yi/JI 

I 

» 

1 1 

, xj 


\ '</{T 

/f/7T 

'>1^1 


\ y. 


(16) 


(17) 


/27 


for the "tripod probe". In both cases, we have: + + ” VT 

i.e. the probe axis is a symmetry axis. 

Building such three-wire probes is surely not much more difficult 
than building X-wire probes. Whether and to what extent they can 
be used to advantage, must also be determined experimentally. But 
it is already obvious that the presence of prongs can greatly invali 
date the analysis in high turbulence, when the instantaneous speed 
direction deviates so much from the probe axis that hot wires lie 
in the wake of the prongs. 
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For non-orthogonal three-wire probes, the system of equations (5) 
generally leads to 2nd, 4th and 8th degree equations for the 
velocity components, so that as a rule, no analytic solution can 
be found. The equation system is simplified for the case of a 
three-wire probe where the wires lie in a plane. 


A special case is discussed here, namely of a plane three-wire 
probe in the form of an isoceles right triangle: 



/28 


We obtain 


^ (w-k^)/2 \ 

4.* 



■1 /I 0 c 

t. 



, 4 »!* 0 / 

^ y 




and, from this: 


' /f 0 0 


“J" ' 




(M 

0-4 4 0 1 

1 * 

, 2 j 

- 


«s 2 

V 1 

1 / 

0 0 0 /!■ j 


* 

i / 

■ 



le-,/ 


in conventional notation 


/29 


+ i ■^4- k*) Oj 


V 


(18.1) 

(18.2) 
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i 

U* t /3 - Cj, . 

(18.3) 

Solving the system 

of equations (18) gives: 



» e, - + £ e, + •/(! 4. + e/ J e/ , 

(19.1) 


- +-■/( 1 ) 

(19.2) 



(19.3) 


The system of equations (18) is invariant to sign reversal of. 

(uj.Uj), u., but not to sign reversal of U 2 or u^ alone. 
Consequently, the probe can be used advantageously in two ways: 


-In sufficiently low turbulence, as for the orthogonal 
2 , axis can be chosen as the main flow direction, fig. 
have : 

U, = +U, =+U.- = ~ u 
I 1 2 3/3 


and 



/ 


U ts 

! 2/3 



probe , the 
5 . We then 


The linear transformation of the base vectors ii> ii is: 



X, / 



4/vr 

^/-/T 0 ^ 


'f 1 

cs 

1 

-'tf-fT + -iZ-Iz I 




1 w/VT 

-yf/fr -'il-fi- 1 


h' 


-In severe turbulence, sign reversal of u-^ is expected. In this 
case, it is better to choose the xj^-axis as the main flow direction, 
so that the probe triangle is perpendicular to the main flow direc- 
tion, fig. 6 . The signs of U 2 and uo can then be determined from 
equation (18.3), if we presume the fluctuations as constant! /30 

in time, if we know the signs at any given time, and if U 2 and U 3 
do not disappear simultaneously. We then obtain:: 


> 

Uj - _ 


(sign u,>o is presumed) 

(sign U 2 , sign U 3 from 18.3 and one 
starting condition) 


Similarly, a plane three-wire probe could be used, whose hot wires 
form an equilateral triangle. The instantaneous speed vectors may 
only lie in the half-space in both cases. 


From a theoretical standpoint, it would be desirable to use 
the least symmetric possible probe, since this provides as much 
information as possible about the flow direction. That the 
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system of equations (5) cannot be resolved analytically in this case, 
need not be an obstacle. If we are satisfied with a numeric solu- 
tion, a fast computer time-signal can give the speed components 
and with the correct sign ( sisjn u^-sign and u^-sign ) . The 

sign Uj^ must be obtained from an initial value and the continuous 
profile of U]^ or by other means (e.g. by an additional hot wire). 

Three-wire probes of different type for an improved measure- 
ment of speeds in highly-turbulent flows were also suggested by 
other authors. Rodi [18] points out that in principle, the instan- 
taneous speed component could be determined from the signals of a 
three-wire probe by means of a digital data analysis. But he 

considers this. method to be very cumbersome and fears strong inter- 
ferences between the three wires. Gaulier [19] uses an orthogonal 
riple-probd'’ by DISA Co., whose wire arrangement resembles the 

orthogonal probe in fig. 3, but with a different prong design. 

Since Gaulier does not use a linearization of hot wire signals, 

he cannot use the advantage of the triple probe over an X-wire 
probe, since he needs a series expansion using small fluctuation 
quantities. In addition, in the analysis, Wrong relations are used 
between the average speed components and averages of the cooling 
rate (eq. (7) in [19]). The same triple probe of Gaulier is used 
by MOffat, Yavuszkurt and Crawford [20] to determine the instantan- 
eous orthogonal speed components in highly-turbulent flows. The 
real time analysis makes use of analog techniques here. The sign 
problem in high turbulence is not discussed. Acrivlellis [21] 
suggests a triple probe as per fig. 6; he also states 'explicitely 
how the square of the speed components :u| = (u^+u])* are calculated 
for linearized signals, but also does not aiscuss the sign reversal 
in higher turbulence. Fabris [22] developed a four-wire probe and 
a method of signal processing and analysis which is supposed to 
permit simultaneous determination of the speed components and 
temperature in turbulent boundary layers with temperature gradients. 
Of the three hot-wires which measure the speed components, two form 
an orthogonal X as in standard X-wire probes, whereas the third 
lies at a 45° angle to the lines of this X, and the fourth — to 
measure temperature fluctuations — is parallel to this perpendicular. 
By using extremely thin platinum— rhodium hot-wires (0.625 jam diame- 
ter), the wires should be "interference-free". If the angle between 
the probe axis and the instantaneous speed is greater than 30°, then 
even for these wires, distortion of signals by the prongs or the 
thicker, silver-coated wire ends cannot be avoided. 


3.3 Two-Dimensional Unsteady Flow 

In isotropic, turbulent flows or in boundary layers, the speed 
fluctuations in the direction of three orthogonal coordinate axes 
are of the same magnitude, so that no component can be neglected. 
But there are flows in which large speed fluctuations occur, not 
because of turbulence, but through other mechanisms and the speed 
vector revolves almost exclusively in a plane. Examples for this 
are the entrainment region of a planar or round free- jet, the flow 
at the outlet of radial impellers in ventilators or compressors, or 
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the flow in the wake of axial-symmetric or two-dimensional bodies. 
These are simultaneously cases in which the instantaneous, local vel- 
ocity vector cannot be restricted to a half-space or quadrant, but 
instantaneous back-flow occurs. Conventional hot-wire measurements 
will have to fail here since the standard formulas for analysis of 
the signals ignore the occurring sign changes of the speed components. 
But it is expected that an accurate hot-wire analysis in two-dimen- 
sions will be much simpler than in three dimensions and that no 
three-wire probes will be needed, but X-wire probes will suffice. 

Let an X-wire (not necessarily orthogonal) lie symmetric to 
the x^-axis according to the assumptions in section 3 : 

•' O <■ 'V / Z • 

fit, wa 0^ ^ j 

and let U3 - 0 . Then for the two signals e^^, B 2 , we have: 



\^2 

<OJ.< 

• • ci>4 0 0 • 

with 


^ -f (k*-- <) 

>*0 ^ 




>0 ^ 


3.^2 

( k* - ) co-i 

«C ' coi ^ ^ 


By 


substitution we 

vA •• 


see that eq. 


= I ( /f 


■ a 


( 21 ) 


has the solutions: 

wc 


( 22 ) 


Now the choice of a, e is basically open: 
a ^ 21 is an orthogonal X-wire probe 

i 

a > 1 is a blunt-angle X-wire probe (related to the opening 
^ angle, with the x^-axis as angle bisector) 

a < j an acute angle X-wire probe L 

g ^ g means that the X lies in the x^X2"plane 

(regardless of the author's reservations ) describes 

an X-wire probe inclined by the angle e to the X]^X2-plane, 

® ^ even though speed components are measured only in this plane. 

®22 is true when cat*a=cos='e thus applies only for blunt or 

right-angle X ( a 2 f ) . In the case where fcot'a > cos“e 
(always true for acute-angle X), then ®ti < ^22 

In general the lines: “2 ~ ° -lines) 

7 = 0 (zero lines of '^■\/2 ) 

tana • cos8u^7 = c (projection of hot wires in 

xiX2-plane 
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do not coincide. 

Example with 


Ml 


< a 


22 




' fow <l> 




. ^ 

Siy* 

cor’^ 0 


( 23 ) 



of • C€»S & 


( 24 ) 


o„ < a 


22 



Key: 1 -projection of the X-wire 


As above, a sign problem occurs: We obtain equations only for 
''1/2^ but not for vj. . In the case of the sign can be deter«i 

mined if it il known at a beginning time and if vi and V£ do 
not disappear simultaneously (or ui and U2 simultaneously). 
Now if we look at 


with 


>/z — vj ' ^ •Sti rr Z 

+■ u/- + 




$ 




( 25 ) 


then we find that this function has extreme values at ^ 4 -^7 
(i.e. '^1= -^2) and no extreme values at ± 4> ± mr, , when <t>/^ . So when 

0 and a are selected so that 4, f , then the following method gives 

the correct sign — at least in principle. Let ''2^° at t=to- If 

is the first zero point following V]^ or V2, then for t^st<t^ 
we also have continuity reasons. If at 
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has no extreme value as a function of the time, then or V£ 
changes sign, otherwise not. A fast computer thus could solve 
equations (22) with the correct sign of Uj^ and U 2 • 

With two-wire probes, in principle the instantaneous values 
of speed of a two-dimensional unsteady flow could be determined. 

The same applies, as shown above with minor restrictions regarding 
sign, for suitable three-wire probes in a 3-dimensional i unsteady/ 
flow. From the time histojry of the speed components, in both cases 
the average speeds and fluctuations of any, high order can be 
determined, even when the fluctuations are very large. The sugges- 
ted probes and the type of signal analysis are however, not in 
practical testing, so it is not known in what cases they are 
practically superior to conventional one and two-wire probes and 
the corresponding signal analysis. As long as one- and two-wire 
probes are used to measure unsteady flows, it is useful to 
consider conventional analysis of hot-wire signals more closely, 
in order to estimate the systematic errors in larger fluctuations. 


4. Conventional Analyses of X-Wire Signals LA— 

In sufficiently small turbulent flows, the average velocity is 
large compared to the velocity fluctuations at nearly all times. In 

this case, from eq. ( 2 ') by a series expansion of the square root, approx- 
imation formulas can be obtained. Such approximations have been 

given by various authors, but many of these series expansions are 
cumbersome. One attempt to derive approximation formulas from a 
consistent series expansion by Bartenwerfer [23], gave in the 
second approximation previously non-standard correction factors 
for the turbulence quantity (i.m= 1,2,3), 

In order to have fewer unknowns, one tries here to set the 
axis intersection so that U 2 = U 3 = 0, or the x^-axis coincides 
with the main flow direction. Finding a criterion for this is the 
first problem; the second probl em t o determine the average velocity 
and the turbulence quantity . In symmetric flows, one will 

naturally take the symmetry into account in the selection of the 
coordinate system: So for axial-symmetric flows without spin, the 
azimuth component of average velocity is zero, likewise the third 
component in two-dimensional flow. 

i • 

If tihe main direction of flow is not known precisely for 
reasons of symmetry, then it must be known at least approximately 
for this method, so that an axis intersection can be 

selected so that the following equation applies to the speed compon- 
ents: 

lu'i ^ Iu,i ; (26) 

Then equation (2') applies in the following form: 
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(27) 


with the abbreviations: 



*-ia 

<a-< 14 


+ 


JilL. 
^/t‘1 u* 



u'/ ’ , uj’ , 

t;/ u/- 




i'J 


+ 


2 ^ iicte + 2 +2-1^ 

■ c5/r^ Oi 




/37 


In c are the first-order terms, and in n* are the second— order 
terms (,,2 not necessarily =■» ) . According to: 

- A -h § X — -^Js ^ 0 (x^) 

the root in (27) can be developed consistently when an is not small 
compared to 1. Normally, these developments (expansions) are given 
only down to the lowest orders for the sake of simplicity. But 
since in the literature sometimes 4th-order terms are found, the 
general expansion is performed to the fourth order as a check. 

With 


Y 2 ^ ' 

we obtain: 


y 


: A -t ^ 4 ^ ( n*“ 7*^ " ^ 

^ . [ ^ + I ■»• i (->s* - f ' 7 - ^ - 


(28) 


^ith 


f 


^ iii ^ 3^ Jff- 
-[T * u, a„,t (4 


^ 2 . i4f + iii Z ■<■ 2 

j ^ \ u/ j. <1 / ^-2^ — ( — )■ + 


/ d« 


wi'* ^/iii ^ - sl-^T) -rh " 

V ( ^44 ' ‘3"- 4 


Since it is already clear that hot wires in various orinetations /3.8 
are needed to determine the speed fluctuations, let us handle only 
the orthogonal X— wire probe below. For such an X— wire, lying 
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symmetric to 

the X]^ axis. 

we thus have 

0 

7T 

1=“2= 4 

and ®2 



( ) 


,, 0-. 

) 

/ 

2 




i 

! -a-ii 


~ 

c 


/ 




i ^42 

( 

V — 


0 .) 

/ 





( ) 







there follows 

from 

(28) : 









»> u ■ 1 

. ^ '5X, TiT 

+ ^ 


Vt 


. ] 




“ <■ 







(29) 


All second-order terms cancel out. Thus, when , by first 

and second approximation, we have: 


Us. ~ ^3 


e, ' e 


•»/* 


(30) 


Thus, the orientation of the X-wire will be varied until iei=e 2 
for all rotation angles 0 . Then, the x.-axis is the main flow 
direction and • 0 ^ more accurately: U 2 /U 1 and U 3 /U 1 are 

then at least on the order of the 3rd order terms, that is, two 
orders of magnitude smaller than the lowest occurring terms. 

One must keep in mind that the criterion (30) was derived for 
a placement of the xn-axis from a series expansion and thus applies 
only for small, perhaps for medium, turbulence intensity. No one 
can say how far this condition is reliable; it can be viewed as a 
calibration of the hot wire. 


4.1 The Second Approximation 

After the alignment of the X— wire into the mean flow direction, 
one can use the series expansion (28) with U 2 =U 3=0 and terms up 
to second order. This seems reasonable since quantities of second / 3 9, 
order, e.g. are to be computed. Since the two signals 

of the X-wire are available simultaneously, the time-averaged 
quadratic quantities- available to calculate the 

correlations and the speed . From (28) follows for one 

X-wire position, the system of equations: 



0 


«A( 3,5 - 

0 


1 0 'i 

A' 

=5 A? 

•9aI 




•2. 






— 2 

Z ^4 

~ 2 


1. 

.9a. 

Ai., aji - .2 


0 
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The second and third equation of the system are precise, 
whereas in the first and fourth equation, terms from third order can 
be neglected. 


By transformation of this [system of equations, we obtain 

2 ( a., Sj 3 - o 

o 0 


0 *" ** ^3^ 






a 

0 

3^ 


Atf^r 

o 




Z ^3^ 




Kr 
T"/ ^ 

U4 Wi 




(32) 



2 . N 


/ 2 . ' 




e^f/a. 





■= 4/m' 

^ Cjf- 




1: t - 1 e»~ 


?( e^' 


J 


s. 


/40 


with et = e. - e.. 

For resolution of velocities, two suitable positions of an 
X-wire probe would suffice. The usual practice dictates selection 
of four orientations, and to be sure, the angles e=o, 7 ,} i, - ff f - 

. 21 . lE since the resolution is simpler and nearly ail correlations 
can be ^obtained from the measured values of one x-wire position. 

a is held constant at f and the corresponding hot-wire signals are 
again designated as e-,. e^, e^, eg, e^, eg 

Thus one obtains a' non^iomogeneous linear system of equations with 
seven unknowns and 16 equations. Not all equations are needed for the 

solution. Equations with the signals e^ to e^ are preferred, so 
we obtain: 
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From equations (33.9) and (33.10), we have: 




2f^-kV ( ' 


(33.1) 

(33.2) 

(33.3) 

(33.4) 

(33.5) 

(33.6) 

(33.7) 

(33.8) 

(33.9) 

(33.10) 

(33.11) 


which is an alternative to equation (33.5) with w hich one can 
make do with the determination of average values e^/g and | e^yg 

Another possibility to calculate is provided by eq. (6.2) with 


U^.U3=0: 


k' 




e« 


+ gf ) 

2. y -» 


and this formula is exactly correct. Equations (6.1) and (6.3) 
show that (33.4) and (33.6) are actually not approximations, but 
precise solutions. 

From eq. (33.7) to (33.10) we also find that in general, the 
average values ®i/ 2 * ® 5 /b' ®?/8 same, and 

this can be confirmed by experiment. Deviations from the average 
values are related to the anisotropy of the turbulence. According 


to (33.7)-(33.9), ^-137 = W ’* 

Apparently, 7^' = if and only if = ‘^ 3 ^ ; at the same time 

is equivalent to " °- 


1 +i<^ 


®i/2 " ®3/^i = 2 ^ 


“ However, e^g* = e, ’ ^ 


Since all averages e^ are measured, one can define anisotropy 
factors: 


and from (17.7), (17.8), (17.2) we obtain: 


UJ 
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(35) 
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with the abbreviation: 

F [ /I + ^ 


/ft- _ , 4-hSk ^ (e«.'-eV;* -j- 


'</Z 


(36) 


For the local turbulence intensities and the normed Reynolds 
shear stresses, there follows: 
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A^-F- 


(37.1) 

(37.2) 

(37.3) 

(37.4) 

(37.5) 


(37.6) 
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4.2 The First Approximation 

In the first approximation, only the linear terms in eq. (28) 
are taken into account and all higher ones are neglected. 


c 



Ut 


. 


] 


(38) 


With the same hot-wire orientations as above, one obtians a linear sys- 
tem of equationsin U. , ul , which can be solved by these variables 
and which then gives'^the^desired time averages. We have: 


fdiU (ea' ej)^ _ -f e/J j- 


ul*- 


1 (&L~ e-7)^ 


'JY- 
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(e; - ] 

J 


(39.1) 

(39.2) /43 

(39.3) 

(39.4) 

(39.5) 
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(39.6) 
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(39.7) 

result 

is obtained if we proceed from 

(38) 

and form 



. There 

results an equation system analogous 
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0 . 0 

O 0 O 
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(40) 




















But a comparison with (31) shows that the relations between 
the speed fluctuations and the hot-wire signals are apparently wrong 
for weak turbulence- By pransf ormation we obtain the system of 
equations analogous to (32): 






o 

0 

0 0 

0 0 


o 

0 




o o o 0 

o o o 0 

0 0 2 <ai, ^3^ 
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+ Ci + 1 ~ 


1 1 1: 


5 e* - Z <5^ C». 


G-f/i 




/ 1 a- 


I (e;* - e'") 

^ C7i~-T /)^ J 


Now the first and second approximations differ only by that 
equation containing the average speed U^, whereas the other equa- 
tions containing only the fluctuation qua(ntities ® . u’ , are iden- 
tical. Thus, via a wrong interim result, we have come to a partly 
correct end-result. But this is not mere chance: If we proceed 
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from a series expansion (28) from 

e = [ ^ * 1 p + O ( J , 

then a linear approximation follows: 

e , 

eT L J / 

e ' •>. e. — S. " (4 ■ ^ * 

- (VSTu,)"- • ^ 

in quadradic approximation: 


e. ' t4 [ "k 


e" « c ^ -f 

f + 1 CV - 

e' « t/SIT 





Apparently, the relations between the quadratic fluctuation 
quantities and ^^e the same in both approximations ,_but not 

e. This means that according to the first approximation, e should 
be independent of e , whereas this does not apply for the second 
approximation. 

Equations (39.1) to (39.7) with k=0 give_ relations which apply 
for low. turbulence and slight tangential cooling; 'these have long been 
known. For greater turbulence and inclusion of tangentxal 
cooling, this first and the second approximation give the following 
corrections : 

-The ”k-correction", expressed by the prefactors .... 

(i + '<^)/(i-i<=‘) in (37). This correction is also provided by 
the first approximation and is easy to see since the coordinate 1 
is assigned the factor 1, and the coordinates 2 and 3 are assigned 
the factor (1+kM/d- Kd . This correction was already given by Cham- 
pagne and Sleicher [24] for some turbulence quantities (eq. 39.2 
and 39.4), if we take into account the identity 

_ 1 

( 1 - j ~ 1 -3k'' +4 k'' 


valid for small k. 

-The "turbulence intensity correction", specified by the term 

. , 3 , je' - y'’ 

~ Z ^ ^ ^ J 

This correction has been given by other authors [25, 26,^27]. 

For turbulence levels up to 0.5, it amounts to about 504. 

-The "anisotropy correction" which is specified by the values for 
the anisotropy coefficients a,u,v . This correction is not provided 
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by the linear approximation. It means that the linear approximation 
is permitted only for the isotropic turbulence (x=u=v=i). 

Figures 7-9 show the correction factors and A* r* as well as 
the corrected average speed for various values of the turbulence 
parameter. For simplicity, the quantities are presented as functions lof 
^ 1/2 “ fixed parameter values of k and > , although 

the anisotropy coefficient ^ surely depends on the turbulence inten- 
sity. For Asi the values coincide According to the relations 

at a plate boundary layer, where applies, when X 2 

is the coordinate perpendicular to the plate, the corrections are 
denoted only for the case A*ii , since in every case the axes x^,X 2 
can be selected so that we have . Estimations from measurements 

show that the factor a generally is almost equal to 1; about o. 95 <A<i /47 

From eq. (37.2) follows : © ^ ^/VT (oak ^ 0 . 2 ) 

when o £ / uiV 0/ ^ o.r 

and we shall limit the discussion to this very region. 

It turns out — and this can be shown analytically — that the 
correction factors fS a^f^ rise continuously with k and ^^ 1/2 , but 
they fall continuously with a. The influence of the tangential 
cooling factor k is small for small turbulence intensity (small 
and weak anisotropy (A=>i) , but grows rapidly with both quantities. 

The influence of turbulence intensity on the correction factors 
is similar in all cases; the curves rise continuously with«^i /2 
and are convex to the bottom. For an anisotropic turbulence,, we 
have f2>a^f* • The anisotropy of the turbulence has a strong in- 

fluence on the correction factors. In particular, x has a great 

effect on F^ . 

5. Analysis of the Quadratic Signals /48 

2 

The attempt to use the quadratic signals e alone for the hot 
wire analysis would eliminate the square root expansion since this apparent- 
ly leads to a systemmatic error in the anal ysi s . Rodi [18] found 
that one can calculate the values up+i7p, from the 

quadratic signals of three suitable hot-wire directions, we 
presume ; but it is not possible to compute the average 

speed U. . In order to determine this, Rodi refers back to the 
series expansion. He supposes that in the case of strong turbulence, 
this method is still better than the conventional one, since 
the series expansion is used only at this point. 

Another attempt to avoid the series expansion comes from Acriv- 
lellis [28] to [31J. He proceeds from the time averages je", 
of differently-oriented hot wires and derives equations for the 
turbulence quantities u'u']; (i,m=i,2.3) , which are supposed to be valid 
for any level of turbulence. This unrestricted validity must be 
doubted since in the course of the analysis, no sign problem occurs; 
and as we know from section 3. 2. 1.1, this is unavoidable for larger 
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turbulences. Actually, Acrivlellis gives a functional relati on 
between the fluctuation signal e*® and the turbulence quantity 
(eq. (21) in [28]), which is not derived and which leads to a 
linear relation e - u (eq. (24) in [28]), regardless of the turbulence 
intensity. From other investigations, it follows that such an 
exact linear relation does not exist ( Bartenwerf er [32]). Thus, 
of the results of Acrivlellis only one is exact and was used in 
the derivation of eq. (21) or (24) — the other results are approxima- 
tions for low turbulence intensity. In a later work [21], Acrivlellis 
also stopped using such a linear relation, in this cas e he did not 
find any exact solution for all sought quantities tiJuX for 

X-wire probes. 


So if an exact analysis with a two-wire probe is not possible, /^9 

then there are yet several possible approximations in the case of 
low turbulence. A very close approximation is to proceed from 
the averages of an X-wire probe, instead of the 

usual thus avoid a square root expansion, but in 

order to get solvable equations, we must neglect terms greater than 
2nd order in t he fl uctuations in the second approximation in the 
terms ! • If u,=u,=o again, or if the criterion (30) 


teTy exists, then one 
to eq. ( 31 ) : 

obtains the following system 
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(42) 


This system of equations is derived in more detail in Appendix D; the 
other transformations and computation steps are presented there also. 
As the final result, we obtain: 

// f I rr 




uLai 

uJ- 


uL 

0/ 


W-k*-/ ^ L ^ J ' 
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• J 
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(43.1) 

(43.2) 

(43.3) 

(43.4) 

( 43 . 5 ) 


/50 


28 



ORIGSNAL FACE IS 

OF POOR QUAUTf 


ui uj 

yf-f- 

4, [ 1 1 i £T t 

(43. 

.6) 
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(43. 
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r y, X i\^ ,) _ <L - 

[ + T^F^'V 2 

(2-f-kVC't+ feV . Ccx~c^)^J- I 


(43.8) 


(43.9) 


If we compare equations (43.1) to (43.9) with equations (37.1) 
to (37.6), (34) and (36), then a strong formal similarity is 
found. But regarding the quadratic signals, note that for turbulent 
signals, etc. applies, thus the appearance of these differences 

in formulas (43.1), (43.2), (43.3), (43.5). 

Also, in this analysis there is a linear approximation. This 
results formally from the second approximation, by setting 
and . The approximation fq='' , but possibly is called 

quasi-linear . 

Measurements with X-wires in turbulent flows could show whether 
the results computed with the quadratic signals differ significantly /5.1 
from the conventional approximations (sec. 4.1 and 4.3). In sec. 6 
these approximations will be included in the numeric testing of 
conventional analysis formulas. 


6. Numeric Check of Some Approximations 

The formulas named in sections 4 and 5 fail when the turbulence 
intensity is high, since in the derivation ''sufficiently'' low tur- 
bulence was presumed. In the measurement, the question is more 
real: How large are the drrors in computed turbulence quantities 
in the various approximations? By experiments it is not 
possible to answer this question since no exact analysis is known 
(e.g. as per sec. 3) which is also reliable in practice. Through 
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a principally different measurement method, e.g. Laser-Doppler 
anemometry, an experimental comparison would be possible; this 
will only be generally accepted as a check of the one measurement 
method if we know which method is more reliable and accurate, 
especially for high turbulence intensity. There is still no 
agreement on this. 


One can suppose that the accuracy of the series expansion 
by small fluctuation quantities is more a mathematic than a physical 
problem and that the used formulas can be tested with simple func- 
tions. These functions do not have to describe real, stochastic 
turbulence signals, they only have to describe correctly the tur- 
bulence intensities in the three coordinate directions, and the ^ 
correlations between the components, including perhaps the bulging 
factor r=u' which plays an important role in other correction 
formulas (Vagt [27]). If these assumptions are correct and if the 
functions are also to be easily integrable, the possibility arises 
of checking formulas used for turbulence quantities by means of 
computer . 


For the periodic step function: 

3 ti', j, 6") - I 


4/s - 


O 4k I' ^ J 


O , 6" * * -T 


(44) 


continued periodically on the whole t' -axis , /53 

with o<s<o<i 



we have (see Appendix E) 


Key: 1-exaraple 
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For the imagined speed fluctuations, we now set: 

y X j V r X , . ^ 

“ u , t(;s - Uj " o. ; 6 (T 6 I ; t random 

In order to use only meaningful parameters ''i« , known experi- 

mental results can be used; here, the works by Elsenaar and Boelsma 
[27], Charnay [28], Dechow [29], and unpublished measurement results 
of Lehmann [30], in order to estimate the value-range of the follow- 
ing turbulence quantities: 




uT- 







(46) 


Although the named papers pertain to entirely different test 
arrangements; they were measured in regions of different turbulence 
intensity and anisotropy, and with hot-wire technology [33-35] or 
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with Laser-Doppler anemometry (LDA) [36] — figures 10 to 16 show 
examples — the following uniform range of parameters results: 


-isotropic turbulence 
-two-dimensional boundary layer 



c, =0 (l*m) ; 
Im 


0.3 S < Cj < 0.6 , 

|C^2l " 0. 3-0.6 . 

^13 ^23 * ’ 


(47.1) 


(47.2) 


-three-dimensional boundary layer 


0.4 < C 2 s < 1.1 . (47.3) 

!c^j = 0 . 0 - 0. 3 . 

~ I‘'i3l ^ 1 ^ 12 ! (with exceptions) 

The x^-axis was always the direction perpendicular to the two- 
dimensional boundary layer. For the turbulence intensity we 
have : 

- 0.001 s < 0.05 for the hot-wire measurements (turbulence 3-24%) 

0.1 < c s 1.2 for the LDA measurements at the round, free- jet 
^ (turbulence 30-110%) 

Accordingly, the following parameters are selected for the 
numeric test: 


-isotropic turbulence 

0 i* * 0. 2 5" ^ 

, 

r,^At 3 ^ 'a “ ^ , 

^ yf/S ^ J • A/C O* ^ 

-two-dimensional boundary layer 


0 « «£ o.2r 

AitA^ - aI/A^ 


o.C 


r, ^ (xf - O. </)/(, A/s -f 9/(a&-3))^ 
Cm A/3 ^ S • Alt 6*w. O.OCrS-^ 

sljit "Ay 

-three-dimensional boundary layer 


r " 5 taw. A0 


A 

c. 


■ 

Cm 


- o.f- , 

- - p. 

/ 

o 

r • 3 ti«.. AO 


(48.1) 


(48.2) 
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C- ® 0. S' 

^ / 
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^«3 ” 3 ( ^ 

r ^ 3 
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Details of the numeric calculation can be found in the programs in 
Appendix F. 


Before presenting the results, three problem areas have to be 
mentioned. The procedure in the numeric test differs from the 
measurement with an X-wire probe since the x^-axis of the former 
can be set as the main flow direction (Uj=Us=o) • With increasing 
turbulence intensity, differences result which point 

up the limitations of the method. The measurement begins — as 
explained above — with the search_for the main flow direction, but 
the criterion for this: is valid only for low turbulence 

Intensity. For greater turbulence, in this manner we do hot find 
the sought flow-line specific coordinate system, rather one with 
average speed components • The turbulence quantities are 

thus not only inaccurately measured with increasing turbulence _ 
Intensity due to the series expansion by the fluctuation quantities, 
but in addition they are described in a "wrong" coordinate system. 
These two sources of error are thus two sides of the same coin, 
namely the approximation nature of the series expansion. 


The step functions used for the test can indeed simulate 
the turbulence intensity and the correlations of the speed compon- 
ents, but the slope depends greatly on the selection 

of parameters s,o (see appendix E). For the used values (s,o) = (i/6, i/3) 
or ( 0 . 0655 , 1/3) it follows (f.!) = (3,o{ or (10,2.6). Whereas for 
the low bulge factor the slope r actually disappears, as 

expected, for the high bulging factor r^iio one obtains a large 
slope for the test functions of 2.6. Now in the series expansions 
the factors for terms of higher order are small compared to 1, 
but appear to be less suitable for such signal analyses, which 
presume a neglect of uneven powers of the fluctuation quantities 
( e. g. Vagt [ 27 ] ) . 

Finally, one could doubt whether the used test functions lead 
to representative error estimations. The supposition also arises 
that the errors in the turbulence quantities are overestimated, 
since in the specified profile of speed, only the extreme values 
are taken, in addition to the averages. The probability distribu- 
tions of the speed components are discrete with only three possible 
values, thus they are not normal distributions. We shall later 
compare our results with those of Tutu and Chevray [16] and Brad- 
bury [17], where normal-distributed speed fluctuations were pre- 
sumed. 


The results of the numeric test are presented in figures 14-24 
as follows: In constant conditions 
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ft * f 

^Jr - - 

various turbulence quantities are plotted as a function of 

the turbulence intensity in the direction of the x^-axis, /u . 
This turbulence intensity varies from 0 to 0.25 or from 0 to 0.05. 
(Greater turbulence intensities than 0.25 were only computed for 
testing, and it was found that in these cases, the computed results 
lead to entirely wrong results). The quadratic quantity 
is plotted as the abscissa, since quadratic quantities Tlpr/uJ are 
to be determined. In the described procedure, the exact quantities 
are linear functions of the intensity and are repre- 

sented in the figure as lines, e.g. 



■ f I 








The error in the various approximations shows up 
from these lines. As relative error, we define 




[( 


Hi 

u 


approx.^ 4xact'' ' ' 


j U 

in 


as a deviation 


exact 


which can be positive or negative. As a second, non linear 
abscissa, the turbulence degree is used, it extends from 


the Interval 0 to 50% or 0 to 22>4 respectively. Now we shall 
discuss individual results. 


Figure 14 shows the relative f luctuationsllp'/u' (i=i,2,3)as a 
function of increasing turbulence intensity o<u^ /u‘<o. 2 S for iso- 
tropic turbulence, as provided by the analysis of the normal 
signal e in first, second and specifically fourth approximation. 

The tangential cooling is neglected here and in the following 
figures (k=0); the bulging, factor is set as r=3. ( The_ inf luence 

of these two factors will be discussed later). For this parameter 
selection results the linear approximation (solid curve) in the 
whole region o<(trp'/uM<o.2S as a good approximation just for this 
quantity , but gives too small values throughout. The second 

approximation (long dashes) gives too large values throughout and 
deviates in the entire interval more greatly from the actual value 
than the first approximation. The special, fourth approximation 
of Vagt [26] (short dashes) proves to be best here; even at a max. 
turbulence intensity of 0.5 it gives a value too small only by 
ca. 1.37o. The relative errors of the first, second and^ fourth 
approximation at (uJ^/u*)«o,i5 are —4.54, +25.54 and —0.04, The 
turbulence intensities and u^/u* prove to be equal in 

the approximations, as must be the case for isotropic turbulence . 

The accuracy of the approximations shows a somewhat different 
trend than for : Of course, the first and fourth approximations 

again give too small values, and the second gives too large values, 
but the second approximation differs from the actual value by less 
than the first one and proves to be the better one here. The 
relative errors at (u^/u* )=o.i5 are —21.14, +3.34 and —1.14. 
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Figure 15 shows the relative fluctuations (’1=1, 2, 3) for 

the same parameter, but this time determined 2 in first or second 
approximation from the "quadratic” signals e^. Apparently, in 
this case the analysis of the quadratic signals leads to much 
greater errors than the analysis of the normal signals e . Of 
course, the second approx, proves to be more accurate than the first 
in the entire range, but for turbulence intensities i^ deviates 

considerably from the true values. The rel. errors at (u^^/u^)=o.i5 
am ount to —57.8% or —18.9% for and —51.1% or — 25.6^ for 

, It will be shown whether the analysis of quadratic 
signals generally lead to greater errors than the analysis of normal 
signals. 


In fig. 16 the rel. fluctuations 1^/u" (1=1, 1,2) are shown 

again; they are obtained from the normal signals and plotted for 
a small range )< 0 . 05 of turbulence intensity (compared to 

fig. 14), thus for turbulence levels 0 to 22.4%. It turns out 
that in this region, the turbulence intensity can be 

defined practically without error by the first or fourth approxi- 
mation, but the turbulence intensities are defined with 

equal accuracy by the second approximation. Figure 17 shows the 
same quantities as fig. 16, but they are computed here by the 
quadratic signals. The second approximation proves to be Useful 
for this region of small to moderate turbulence, but stays clearly 
behind the analysis of normal signals. 


Figure 18 shows the same quantities as fig. 14, but now for the 
case of non- isotropic turbulence. The selection of parameters 
0 , 3 . “|”-r/17Ty=o.5~ tT|TJ/(#p’.VS7)= -o.4. ■np7J/(/r|7.V^) , 0 

should simulate a two-dimensional boundary layer. In comparison 
to fig. 14, all appro^mations of the "normal" analysis for tur-_ 
bulence intensities o<u^^/ u^SD‘ 10 prove to be just about as good as in 
the case of isotropic turbulence. Above ilp'/u^=o.io to 0.15, in all 
approximations a clear bending off from the exact l ine ar profile 
occurs, which even leads to negative gradients for 
to very large errors for all turbulence quantities. For o<u^ /u* 
the second and fourth approximations are more accurate than the 
first, } however, not- to the same extent as for i sotropi c turbulence. 
Figure 19 again shows the turbulence intensities ujVuMi=i . 2 , 3 ) 
computed with the linear and second approximation of the quadratic 
signals. Whereas the first approx, proves to be definitely in- 
accurate, the second one lies surprisingly close to the true vAlues 
in contrast to the case of isotropic turbulence. A plausible 
explanation cannot be given for this. 

Figure 20 shows computed approximation values for the sole 


shear stress quantity 




/U^ 


which does not disappear in the 


two-dimensional boundary layer. In the upper part of the figure 
we see that the normal signal leads to good approximation values, 
as long as L(j*/u^<o.io . As in fig. 18 for the turbulence intensities 
, for values above (up/u*)»: 0 .iD to 0.15 there results a clear ^ 
bending of the approx, curves from the exact lines. The analysis 
of the quadratic signals — middle part of the figure — only gives good 
results here for the first approximation. The lower part of fig. 18 
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shows that the striking deterioration of various approximations 
above u’p’/u*>o.i5 is related to the approximation character of the 

calibration criterion for finding the average flow direction; see 
section 4. In contrast to measurements, here we have_ set U2=U3=0. 
From a certain turbulence intensity, a difference of results. 

Thus it can be concluded that above this certain value of turbulence 
intensity, the specification of the flow— related coordinate system 
i(u =u>o, U 2 =li,=o) is problematic. 

From eq. (6.1) we can derive the relation between and 
the difference of averages • For U 2 =Uj=o follows form (6.1): 






and again: 


TTTi 2 { r __ 

“ * (. 2 ) { e, + e* i 




e, + c, J 


If we divide by the average speed U, we have: 
\ O’- fexact L ^ o* Misipp. e, + Ci J 


/6l 


i ^ ) 

(F from eq. (36)).j \ o*- .,app. 

So if e^>e 2 , then the 2nd approx, gives too large values, for 
<62 too small. Figure 20 confirms this and also the quantitative 
relation if we note that here 'u^u^<o . 

In the case of a three-dimensional boundary layer, all correla- 
tions tTfC''" (J 7 m) are different from zer o_^„., In,^e follov^ng figures 

the param'eters uT^/Lrp‘=o.5, irp'/u|»=D.7, u*uJ/Vu|*' u^= 

and ^jjuj/-fu^u^-. - 0.13 were selected. Since the analysis of the quad- 
ratic signals leads to worse results than the analysis of normal 
signals almost everywhere, only results of the latter are shown. 
Figure 21 gives the turbulence Intensities and m first 

to fourth approximations. The various approximations do not 
differ significantly here. The curves are similar to the case ot 
2-dim. boundary layer (see fig. 18), especially the same bending 
of the curves from the exact linear run in the range of turbulence 
intensity(i;F/u")'*’0.10-0.15 is found. The errors above 0.15 are about 
the same magnitude and below 0.10, they are slightly larger than for 
the 2-dim. boundary layer. The same thing applies for the shear 
stress quantity fig- 22, upper section. The approximations 

differ only a little from each other, but exhibit very large errors 
above 0.10. Above rrp'/iF=o. 175 all approximations for 

this parameter selection give a wrong sign. For all turbulence 
intensities, the error in is greater than for the 2-dim. 

boundary layer. The lower part of fig* 22 shows u^u^/u and the 
first approximation gives improbably accurate values, even fot 
extremely high turbulence intensity of 0.25 (= 504 turbulence level). 
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But this result cannot be generalized, since it is a res ult of 
the special selection of parameters. The profile of U 2 U 3 /U 
over the turbulence intensity u|^/u^ , fig. 23 top, is again^ 

"normal”: For turbulence intensities below 0.10 good approxima- 
tions, above 0.10 large errors increasing quickly with the turbulence 
intensity. The lower part of the figure again shows 
representative for the other dif ferences j )/( 5 . ^5 i» 3 , 5 , 7 . 

whose profile is quite similar. As in the case of a two-dimensional 
boundary layer, we see a sudden rise in these differences in the 
same region I]tt/u>»o.io-q. is in which the frequent bending of the 
graph (mentioned above) of different approximations from the 
exact profile occurs. 

If we summarize previous results, then we find: 

—All examined analyses of hot-wire signals based on series expan- 
sion by small fluctuations, lead to systematic errors in all 
turbulence quantities FJlPJ/u’ (i,m=i,2,3) for a high turbulence inten- 
sity 


-The analysis of the "squared" signals usually leads to worse 
results than the analysis of "normal" signals and can thus be 
discarded. 

The systematic error is linked with the approxij^mapion character 
of the criterion for finding the direction of the average speed 
and to specify a "flow— line specific" coordinate system. 

—Xhe turbulence intensities ujr/u* (i=i, 2 , 3 )are a s a rule computed with 
a smaller error than the shear stress terms if these 
are different from zero. But this was already expected, since 
the latter terms are generally almost an order of magnitude 
smaller than the former. 

-The different approximations are loaded with errors of comparable 
magnitude for high turbulence. In particular, it is not true that 
the higher approximation is generally the better one. 

-One can distinguish three regions of ascending turbulence inten- 
sity (naturally somewhat arbitrarily within the limits) in which 
the approximations are good, useable and useless: 

Region irp/u' 

— ~1 ■(T-'S' . 'Us GTTI3 BTiTTS 

II 0.05-0.10 22-32 0-50 O-Oher Key: uber=more 

III > 0.10 > 32 O-iibsr 100 J 100 than 


—For regions II and III? or regions of high and very high turbu 
lence intensity, the deviations of approximation values from the 
actual values can become very great and depend on the turbulence 
structure without any discernable trend. In region I, the devia- 
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tions increase slowly with the turbulence intensity by a minimum 
quadratic amount. 

Since in many cases cf turbulent flows, the turbulence intensity 
is less than 20% and since in this range the systematic error 
increases by the series expansion by the small fluctuation quanti- 
ties in a monotonous manner with the turbulence intensity, then 
by consideration of as many results as possible we could check 
whether the interaction of turbulence structure and tangential 
cooling affect the error in a dlscernable manner. For this purpose, 
in tables 1 to 3 the errors in the turbulence intensity iiy57u*=o.05 
(turbulence degree 22.4%) are presented for a large number of 
parameter values of k,r, up'/u^, up /up and an example 

fig. 24 also shows the influence of tangential cooling and of the 
bulging (convexity) factor r on up/u^ for region I. In this 
example, the influence of tangential cooling is low, but that of 
the bulging factor is great. 

Tables 1 to 3 point up trends, but no Invariant rules. In 
general, the studied systematic error increases with increasing 
anisotropy of the turbulence, and also with increasing bulging 
factor r. The tangential cooling changes the error only a little, 
but it can be larger or smaller. The result presented above that 
the turbulence intensities are determined more accurately than the 
shear stress terms is also confirmed. If the latter disappear, 
then the approximations give the value zero, even at high turbulence. 
But for finite values, the error can apparently be up to 95%, and 
the bulging factor does not have to be large. In the example with 
this large^^error (table 3, 3rd and 4th line), a 3-dimensional 
boundary layer is simulated with 

<= - 0.10 , 

= - 0.10 , 

= 0.27 . 

The large errors occur in the quantities u|uj/u^ and , which 

are determined one order of magnitude too small. At even greater 
turbulence intensity, a wrong sign even occurs and even greater 
errors. 

For small bulging factors {F23) in the case of isotropic tur- 
bulence at the point Lp/u“=D.05 , the maximum error of turbulence 

intensities 'p/u*' for the first, second and fourth approximations 
is 12, 10 and 67o; for very high bulging factor (r»io) 20%., 31%, and 15%,. 
For anisotropic turbulence (two- and three-dimensional boundary 
layer) the corresponding max. errors are 14%, 16%, and 14% (r=3) 
or 32%, 29% and 31% (r»io) . (The large errors ^33 given in the table 
for the fourth approximation were not taken into account because 
they arise from the incorrect assumption up=up ). For the 2-D 
boundary layer, the max. error at low ‘^and 61% at large 

T. In the 3-D boundary layer finally, the named max. error of 
around occurs even for small bulging (convexity) factor. 
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As a final result we find: 


-For turbule nce intensities between 0 and 20%, the turbulence 
quantities u]u^/u^ are determined with an error which depends on 
the turbulence structure with no discernable trend. Thus no 
method can be given for correcting this error. For the shear 
stress term, this error can become so large that its determination 
at turbulence intensities above 20 %> must be placed in question. 

Tutu and Chevray [16] showed in 1975 that in the determinati_on 
of the average speed U and of the turbulent fluctuations 
in 2-D boundary or shear layers via X-wire probes in high turbulence, 
very large errors can occur. These errors, it became clear, are 
based not only on the termination of the Taylor series expansion by 
the small fluctuations, but also on the following rectification 
effect of the X-wire probe: When neglecting the tangential cooling 
(k=s 0 ) and the speed component U 3 , the following expression from 
(21) applies for the orthogonal X-wire: 

e*, = i(ui ± Ui ) 

j /2 

Ot 1 I 1 

ewj = ~ |ui ± ujI . 

/2 

But in the signal analysis, if we proceed from the relation 

ei/2 = (uj i U2) 

I ■ , ' 

which proceeds from the linear approximation, then we preslime that 
due to ej /2 >0 we also have: 

- ^ < atctan(uj/ui ) < ^ . 

This condition is not fulfilled in higher turbulence. 

Now Tutu and Chevray calculated the errors in /66 

u = Ui, 

assuming normal-distributed probability distributions of the speed 
components as functions of for different parameters 

, k ^nd presented them in a table (table 1 in [16]). 

By using these values, one can set up the following comparison table: 



k = 0, 

Ci2 * Cj: 

= C J = c 

.64 


l± 


A j j C ? ] 


Aa 2 r ^ 3 


Aj al %1 

IJ^ 1 

1 

nach [1 6j Autor > 

<L 

nach f 1 B) Autor t 

nach f1 6] Autor , 




1 .NSh. 


1 ..Mah. 


1 .iMah. 

0.01 

10 

-1.G 

-0.5 

-2.6 

-1 .0 

-2.6 

-0.96 

0.04 

20 

-7.B 

-2.4 

-10.7 

-4.B 

-11.5 

-4.3 

D.09 

30 

-18. B 

-4.7 

-24,8 

-B.B 

-29.5 

-8.6 

0.16 

40 

-32.9 

..33.7 

-41.2 

-42.5 

-49.3 

-29.9 

0.25 

50 

-46,1 

-57.6 

-55.0 

-63.0 

-64.2 

-48.1 


Key: 


a-from [ 16 ] , 
1 st approx. 
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Among these parameters the errors are apparently comparable. 

A precise agreement cannot be expected, since in both investiga- 
tions it was found that the errors depend a great deal on the tur- 
bulence structure. Even the fast increase in errors for turbulence 
degrees of 20-307o shows agreement. 


A corresponding error estimation (to the one of Tutu and 
Chevray) was conducted by Bradbury [17] for the one-wire probe. 
He also proceeded from normal-distributed fluctuations, but 
plots 


U ^/U j and 

exact measured 


uJ 


T~ 

exact 



T~ 

measured 


as a function of 


767 


u 


7T- 
1 

measured 


/u 


measured 


for various parameters 


u; 


1 — 
exact 



exact 
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In this comparison too, it turns out that the errors are of 
the same magnitude, even though one-wire probes are compared with 
two-wire probes; the errors stay within a reasonable range only 
below a turbulence degree of 30-40%, that is, they stay below 40%.. 

The error estimations of Tutu and Chevray and of Bradbury 
relate only to the isotropic turbulence or to the two-dimensional 
boundary layer. For these cases, errors resulted in the present 
investigation which are on average smaller than in the case of 
the 3-D boundary layer and which are otherwise of the same magni- 
tude as in [16, 17]. This result supports the contention that 
representative errors are estimated here. 

7. The "Undirected” X-Wire 

It is naturally troublesome, but generally unavoidable, to 
search for the main flow direction for each measured point, to 
compute the turbulence quantities with respect to a flow-line 
related coordinate system, and finally to transform it into a 
fixed coordinate system. If the main flow direction changes only 
a little, e.g. in a round free- jet, in which it is everywhere 
nearly parallel to the jet axis, the analysis of sec. 4 can be 
generalized to a fixed coordinate system. 

Now when the main flow direction is nearly constant, then the 
x.-axis can be fixed (in the round free- jet, e.g. as jet axis) so 
that : 

» lu^l, IU3I 


and for sufficiently low turbulence: 

Thus, one can proceed again from the expansion (28), but now an 
alignment in the local main flow direction can be omitted. 

As in section 4, there are again approximations of varying order 
depending on how far the series expansion is carried. Since the 
first approximation has proven to be good for low turbulence, as 
a rule the first approximation will also suffice here when j 

and I are sufficiently small. At U2=0 the values 

lu^/u^ I = 0.1 , 0.2 , 0.3 

correspond I. to the deviations at 5.7°, 11.3°, 16.7° of the main 
flow direction from the x.-axis. For greater angles, the second 
expansion should be taken"^ into consideration; this is presented 
in appendix G. 


In close approximation, for an X-wire lying symmetric to the 
x^-axis, we have: 





A + 


( J 4 , i:! 1 + ^ 4. fdi ^ 1 
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For an orthogonal X-wire in the usual positions 
e=o, 7 T; n/ 2 . -n/ 2 : ii/A. - 3 TT/ 4 ; -ti/4, 3tt/ 4 with the signals e^ to eg, it follows 

PG-T^ r 1 1 L 3 , 


^‘iit 

«3/y 




[ ( (J4 +Uj) + 

.1, .. I, ~ 


From this we obtain: 




(e* -LeZll 


H,u, "■ 


- 


14 


/( + (f ^ 


2 

/ yi* 

yt-f- *r*- 

\ 

Z 

t y1 + k^ \2 . 

yf+K^ 

( ^-Ar* > 

2 

i ■ 

yf-tk^ 


Z 

t yU-k*' \2. f 

yi+H^ 

V /f- J [ 

2 

/ yl-hk'- \ * 

yl + k’- 


1 ^ 



/IH-K’- 


' ■/ 


Z 

yi-h 

yl~ Ic*- 
xf + fe*- 


■i" 

- eJ^y 


( ej-es')’ 

4- 


(ep - eV>^ 




Cl. - fc„ 
e* 


(49.1) 

(49.2) 

(49.3) 

(49.4) 

(49.5) 

(49.6) 

(49.7) 

(49.8) 

(49.9) 


y>- K» 


The equations (49.1) to (49.6) are identical with (33.1) to 
(33.6); so in the signal evaluation of the fluctuations, nothing 
is changed. The components of the middle flow are given by (49.7) 
to (49.9), where necessarily the following expression must apply: 

;l ~ /l^» +cl*( ^ I Cv - 5,1/ I eV + e, / w 

Proceeding from equations (6.1) to (6.4) and the linearized 
equations above, by neglecting the fluctuations, the error in the 
velocity components U^, U 2 , Ug can be accurately determined due to 
the linearization. We have: 


with 



1 f 

S ’'***«* 1 


1 

(U^ jlin. 

7^ 1 
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for “ 0 . 
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For the example lIj/u, = o. u^/u, = 0 . 1 , 0 . 2 , 0,3 we obtain the following 
number values: 
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u./u. 

angle 

devia- 

tion 

1 

X2 

X 

0.1 

5.7° 

0.900 

1 .001 

0.2 

11.3° 

0.927 

1 .039 

0.3 

16.7° 

0.652 

1 .083 


The large component is thus determined as too large in the 
linear approximation, the small components U 2 ,U 3 are too small 
and have about twice the deviation as the component U^. For tur- 
bulent flows, these errors can be even greater. 

Practically speaking, it is also of interest to know how much 
the turbulence quantities in this coordinate system differ 

from the corresponding correlations arbitrary flow-line 

related coordinate system (V 2 =V 3 = 0 ). The question is usually how 
far can a transformation of measured quantities be omitted; up to 
what deviation of coordinate systems is this possible, or at what 
ratios U 2 /U 1 , U 3 /U 1 . This question will be answered in the following 
section. 


8 . Transformation of Fluctuation Quantities 
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Two ortho-normal, right-hand-dri ented base vector systems 
(coordinate systems) and ( 1 = 1 , 2 , 3) can be converted one into 
the other by rotation: 


^ ( Amj) 




'Mi 




»i 1 




(50.1) 

(50.2) 

(50.3) 


The matrix elements 
of the corresponding base 
angle: . 


of matrix A are the scalar products 
vectors, or the cosine of the enclosed 


? A velocity vector u will have the components u^^ or v^ with respect 
to these bases : ^ ^ 

u " if, =• £ V V, 


We then have : 


or 

Vj « y ■ pi - 
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or briefly: 
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to 


Since A is a linear, time-independent operator, the same 

relation applies both for the average velocity components U^, as 
for the fluctuations u|, v| : 


U •= A ^ V 
U* A. • V * 

use “■ 


(51.1) 

(51.2) 


Consequently, we have; 


V, 

2 S- 



(52.1) 
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i/« u„ ^ 

(52.2) 

(52.3) 
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Three examples: 
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(53) 


By means of eq. (53) we now also find how much the turbulence 
quantities differ from the for small angles y , or when 
the transformation can still be omitted. 


Equation (53) is of interest only in a non-isotropic case. 
Let us consider the 2D boundary layer, and take only the shear 
stresses: 
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SJ 


lAI 


cos 2 y 
cos 2 y 



With the values gained from experience: 


^ fsc o. 5” 


there follows: 


^ ^ coj- + -^s.-p^Zy 

h>7 

Examples: 


Hi 

Ui 


zm . 1 

u'ui 

0.01 

o.s?' 

-0.014 

0,03 

1.72“ 

-0,044 

0.05 

2.86“ 

-0.075 

0.1 

5.7“ 

-0.158 

0.2 

11.3“ 

-0.346 

0.3 • 

16.7“ 

-0,550 


Apparently, we will only be able to omit the transformation 
for very small angles since a difference of 4.4^ shows up even 
for an angle of 1.7°, and it increases more than linearly for larger 
angles. And it must not be forgotten that due to the linearization, 
the angle y is underestimated (see sec. 7). 

-For the orthogonal probe from fig. 3, we have: 
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I ylt-fW ® 2(-Ti 

\ ^ I'Tf /f/Vj 



Thus : 


I yit-n ^/«/T \ 

A = 0 -'Urz yil-fi I 

\ 2/VT -ml ' 


where x-, (1=1, 2, 3) is the orthogonal coordinate system given by the 
hot-wire directions, whereas % coordinate system defined by 
the probe axis and two directions orthogonal to it (as indicated in 
fig. 3). Thus we have: 
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-For the probe design as per fig. 4, we have 


A 


A 


' Al fi 

~2/rr 

0 

At-n 

Ai-a 

-A/fi 


A /ft 

ytl-Ji. 

1 At-IJ 

-A/-/3 

Aj-fi 

~2/-a 

A/AZ 


\ ^ 

-A//2 

At-rz 




Hi‘^1 




(54.1) 


and 
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(54.2) 


Theoretically, the transformation of the fluctuation quantities 
is quite simple, but practically the following problems come up: 

-Even if only one fluctuation computed, in general 

the entire set "" 


u u must be known. 

I m 


increases. 


-Due to the linear transformation, the _ error in . - - 

since the measuring errors of the individual, included 

are additive. 

9. Correlations of Frequency Components 

The frequency spectrum of the speed fluctuations of a turbulent 
flow is broad; in addition, it can contain discrete frequency com- 
ponents, e.g. for periodic, external disturbances^ or in rotating 
flow machines. Thus we will be interested not only in the correla 
tions of the fluctuations themselves, but also in the correia 

tions u' u ~ *" ' of individual frequency components. As' long as linear 

lu) inu) 

relations prevail between the speed fluctuations and certain hot-wire 
signals — e.g. in the exact analysis (sec. 3) or in the 1st approxi- 
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mation (sec. 4.2) — the equations derived there apply equally for 
individual frequency components since eq. (14) or (38) can be 
Fourier-transformed entirely. But note that for the 3-wire probe, 
the signals e. are not tp be filtered, but the signals ^.! 

For the second approximation it was found that the equations 

(33.1) to (33.6) result for the fluctuation quantities and 

from the linear expansion (38), but not equations (33.7) to (33.10), 
which contain the average speed Consequently, the equations 

(37>.l) to (37.6) of the second approximation also apply for fre- 
quency components if the unfiltered fluctuation quantities are 
used for the average values and for the correction factors. 


10. The One-Wire Probe 


A simple, single hot-wire is used for speed measurements more 
frequently than an X-wire, e.g. when we only want the average 
speed U and the turbulence intensity TTpVu^ or when we only want 
a qualitative picture of the flow. But what can be said about 
the accuracy of these measurements? 


The main direction of flow is found by the criterion (30), 
but not in a manner so that the probe rotates continually and 
the signal must be symmetric to the angle a=f . Next, the main 
direction of flow is chosen as x^-axis (U=U^) and the probe is 
used in the position : 


(,u + ]■ + 


(55) 


Although in this case we have only four variables for which to 
solve — (u+u^h ui,^ — this solution is not possible no 

matter how many values e we use. We obtian only . one exact equa- 
tion: 

2 ( co£ ^ (for any ® ) (56) 


This quantity is of interest in measurements with one-wire 
probes in rare cases only. To compute the other quantities u. , 
approximations are needed agian. A consistent series expansion 
out to the second order (as in sec. 4.1) gives: 
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(57.1) 

(57.2) 

(57.3) 
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The equations are not sufficient for computing li, and 
From (57.2) and (57.3) we know only: 


In order to determine the average velocity U, we set: 




u 


, 3 - 

a/3 


(^zt3 ' 


(59) 


with values gained from experience C 2 73 of 0.5 (boundary layer) 
to 1.0 (isotropic turbulence). Then we have: 

^ e;* • [ ^ eiVe,^ j - 

SB • [ /i i ^3 "f" 

uT*- ^ — « 

ijx yf — C Cz'f' C3) e.^ I Cj. 

_1 (61) 
,, 

” yj ^ Cj ■/ 


A linear approximation gives: 


and 



(62) 

(63) 


According to the results of section 6 for the X-wire probe, it 
should be expected that even for the one-wire probe, the second 
approximation will not give much more accurate values than the 
first approx, and that the deviations of these approx, values from 
the actual ones will lie on the same order of magnitude. 


11. Slowly Rotating Hot-Wire Probe 
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Instead of a one-wire probe set into different positions in 
sequence, it can also be rotated about its axis so that the hot 
wire is always at the same place. The rotation should be slow 
enough that during the time averaging of the signal, the orientation 
of the wire can be viewed as constant. 


From eq. (4) it follows: 

^ ' 2 nr 

c= [ -/] ( 

+ £ [ 2 Ch’- '"*> ^ \ ^ 

, / r 2 coj «• •S''*’”' 3 ^ ►> 0 
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specifically for: 


o( ^ T/2 


e* 


o< » v/v 


(normal hot-wire probe) 
(slant -wire probe) 




3*k^ 


34A’ 


+ + -9- «, 


+ •cojg' + (fe*--^) U,l^3 ' + 

^ ( c<^ _ • COJ 2© -t- t<a.<43 • -Hin 2 & . 

Through an analysis of the fun ctional dependence ‘P'i.e) , at 
most the correlations u j' u* , uiu$, ujus can be determined, and the 
difference i“f - uT and a linear combination of but not the 

individual terms iT» “T separately. This result thus 

corresponds exactly to the result of section 3, equations (6.1) to 
(6.4). 

Let us consider the case of low turbulence intensity and 
presume that the main direction of flow is known, thus: 


(4 


0> 


and 


lu: I « U, 


Now we can distinguish two essentially different cases, namely 
that the rotation axis (= probe axis) is aligned parallel or 
perpendicular to the main flow. 

If the rotation axis is parallel to the main flow direction, 
then from equations (31) or (40), we have the following expression 
both for the second and first approximation: 
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f a 




( - s:") 


« ‘-0*' 


^ 4- 




+ 2 Uj^ t(j 


. T ^ r ‘ ‘ 

4 - [ Ck*'“'1) cos c< ■ smci J Kj ■ sin 




$ 4 - 


• Sin & + 


4 - 2 + C.OS 1 . sin’ll- <'<4 *^ 2 . • 

^2. [^ + iK’--3)c«/.4 3 CCK’-3; C4>f.<s.v,«] 

+ 2 [ cw« <^Z <^3 • Cor <9- ^ 


(65) 
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with the abbreviations: 



2 t2 . 

C .01 X J 



(66.1) 


= 1 

coix 

lT^ 

/ 


(66.2) 

Wj 

^ [ 

^ Z 

toi X JimXJ 

■ ^3 / 


(66.3) 

Wiy 




(66.4) 

7 

Iff 2. t / 



(66.5) 

wu 

- 2[ 


Uj 0,3' 


(66.6) 


a,. is 

expansion: 


thus a periodic function of e with 
2. 

Q e'^ t ® ^ *“ 4 ^ C.OS \> B + &tr ■ 


a 




finite 

s>0) . 


Fourier 


A very careful determination of the Fourier coefficients 
gives the values: i > 

° UU + ^{U2-*W3) = “o 

U)2 -U J = 2 «2 

bf ii - Ct 1 

U5 « B.l 

U 6 = 2 62 

A normal one— wire probe ( 0 = 11 / 2 ) would only permit determination 
of "uP" ; a slant-wire probe will ,give in principle 

-the correlations ujir'' , TTHiy , uTuT . 

12 2 9 


-the difference 




-the linear combination 


[ /U- ik^-yi) CoS^«t] ^ 4 




but not the individual turbulence quantities Thus the 

result resembles the general case, equation (64), quite closely. 
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Things look better if in the case of low turbulence intensity 
the rotation axis stands perpendicular, or nearly so, to the main 
flow direction. In order to retain 9 as rotation angle (for «=const), 
the polar angles are now selected with respect to the X 2 ~axis 

(= rotation axis) and 1^1 1 . |u,+u^|, |U34uJi« u, are presumed. In linear 
approximation, we obtain: 


e 



4- U j 1 


(67) 
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and 322 ”^^7 be small. If we choose for instance, then 

> 1/2 for all e. 

For the average signal, after brief intermediate steps, we have: 
FaTT e » 1.4-f- J • U, + [ Uj coj <? + 

( 68 ) 

+ [ ^"i" l ' \J^ Cas2& + [ Ji'n^c<] L>2 Sin2& . 


From this, U. , U 2 , U 3 can be determined and the prerequisite 
|u*l,|Usl « Ui must be confirmed. 


For the fluctuation quantities, from (67) we have 


•*Z’l 




a ^ X. ^ iaJ'^ 




+- 2 ax, uJmx f 2 K u, + 


+ 2 ^xx a^g 


(69) 


After lengthy calculation it turns out that agz as a function 

of « now has a Fourier expansion: 

with the coefficients: 

«, [ /t + \ I + (70.1) 



■= f Sin 2.d ( \ )3 ^ 

(70.2) 


- [ cojrf ] . Ui Uj ^ 

(70.3) 

0(. •2_ 

•= [ (k^~") S'n^^< ( ^ + ^ + 

(70.4) 


+ [ "I ^ 



•ot Si > 7^*)^ yf -f~ — ^ X/W^of ) J y 

(70.5) 

<<3 

f ^ X#'« 2.*y 3 ' / 

(70.6) 

h 

■'*' [ ^ JT/b^ •< ■ J<« 2f< ] ■ Ki L<3^ ^ 

(70.7) 


- [ f ] -(^ - u7*) ^ 

(70.8) 


- [ ^ J U,V ■ 

f 70. Q) 
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The quantities and FJTT]’ result from eq. (70.5) or (70.9), 

(70.2) or (70.6), (70.3) or (70.7) directly, when 0 < a < n/ 2 . . 

The equations (70.1, 70.4 and 70.8) give a linear system of equations, 
in which ipnr ■jprr, TTTT are always solvable for k « 1 and 0 < a < n/z . 

1*23 
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De Grande and Kool [37] came to the same result for the case of a 
non-linearized, slanting one-wire probe and small fluctuations. 

When “ = "/2 , or for a normal probe, from equations (70.1) to 
(70.9) only the quantities uTuX (from (70.5) or (70.9)) and ^ /85 

can be determined. This result too, was found earlier for a non- 
linear probe and small fluctuations, by Fujita and Kovasznay [38]. 

In addition, it became clear here that even for slowly rotating 
probes, a precise signal analysis is not’ possible , rather a series 
expansion by the small fluctuation quantities is unavoidable. Thus 
large errors result with this method too, for large turbulence 
intensity. 
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12. Summary /86 

The measurement of flow speeds and in turbulent flows, of the 
speed fluctuations by means of hot-wire probes and anemometers, 
is a widely used and indespensible technology, although today 
various contactless (probe-free) laser measuring methods are avail- 
able. In the determination of speed fluctuations and correlations 
of various speed components, the hot-wire measuring techniques 
with multi-wire probes are by no means inferior to the laser methods, 
but they are less complicated. 

Through decades of development, the hot-wire anemometry has 
been constantly improved and its useage expanded. But like all 
measuring techniques, it has its limitations which have to be 
known for a realistic error estimation in measured values. 

Besides possible errors due to the electronic control of the ane- 
mometer and due to interactions between the probe and the flow, 
the usual methods of signal analysis for one- and two-wire probes 
contain a systematic error, since the functional relation between _ 
the hot wire signals and the components of the velocity 3-te approxiir.a— 
ted by a series expansion by the velocity fluctuations. Naturally 
it must be assumed that the components of the velocity fluctuations 
are small compared to the average velocity. This error is as a rule 
even greater, the higher the level of turbulence. Above a certain 
turbulence level, the use of conventional hot-wire anemometers is 
no longer meaningful — regardless of other possible errors which 
can also increase with the turbulence intensity. 

A precise consideration of linearized hot-wire signals leads 
to the conclusion that hot-wire probes with at least three inde- 
pendent hot wires will permit an accurate signal analysis in the 
sense that the time history profile of the speed can be determined 
by direction and magnitude from the time history of the hot-wire 
signals. From three orthogonal components u^ ( t ) , 1=1, 2, 3, for /°7 

example, both the components of average speed and random averages 
of the fluctuations can be determined, including the turbulence 
intensities up"* i = i. 2 , 3, and the correlations Lpjj;, i^m. The 

direction vectors determining the hot-wire orientations need not 
span the 3-dimensional space, but can simply point to 3 different 
directions in a plane. 

For the obvious case of three perpendicular hot-wires, the 
analytic relation between the signals and the corresponding orthog- 
onal speed components — a quadratic form in the speed components — 
is very simple, but due to the high level of symmetry of such an 
arrangement, information is lost. Since the quadratic form is 
diagonal in this case, we obtain only the square of the speed com- 
ponents Uj^(t), 1 = 1 , 2, 3 and in case the sign of the speed component 
changes over time, the particular, correct sign must be found by 
other means, e.g- by other hot-wires. If occurring sign changes 
are not taken into account, then too high averages and too small 
turbulence intensities are measured. 
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For the other extreme case of an asymmetric three-wire probe, 
the resolution of the quadratic shape into the spedd components 
is unique (except for one sign), but it is not analytically possible, 
so that one then has to rely on numeric methods. For the case of 
a "smooth” 3 -wire probe, i.e. when the direction vectors of the 
three wires lie in a plane, different examples can be found for 
which the resolution of the square shape is simple and the symmetry 
of the probe is less than for the orthogonal probe. Depending on 
the type of turbulent flow, it may be useful to set this plane 
perpendicular to the main flow direction or at a specified acute 
angle. 

The directions of the hot-wires does not specify their positions, 
so that even different wire arrangements are possible. In sections 
3.2.1 and 3.2.2 examples are presented. Experiments will have to 
show which probe design is suitable for which flow, since inter- 
actions between 'the wires and reactions of prongs and shaft on the 
wire flow cannot be described theoretically. 

In two-dimensional flows in which a speed component can be 
neglected (not only its fluctuations), e.g. in the wake of axial- 
symmetric or 2 -dimensional bodies, the two remaining speed compon- 
ents can be accurately determined from the signals of conventional 
X-wire probes. Again, to avoid ambiguity it is useful not to have 
the probe at max. symmetry. One can use either non-orthogonal X- 
wire probes or the flat, conventional orthogonal X-wire probes 
screwed slightly out of the speed plane. In this manner, speeds . 
and their fluctuations can be measured in regions in which instan- 
taneous reactions and extremely high local turbulence intensities 
occur. 

As long as hot-wire anemometers with three-wire probes are 
not adequately tested and reliable, in the future turbulent flows 
will continue to be measured with standard orthogonal X-wire probes, 
even in high levels of turbulence. The different conventional 
analyses of such hot-wire signals are formally very similar. First, 
an orthogonal coordinate system is defined by means of the X— wire 
probe, such that Uo=U 3 = 0 . (This "flow-line related" coordinate 
system will generally be different at each point of the flow field) . 
Next, from the hot-wir^ signals , the average speed U=U^ and the 
turbulence quantities are determined. In both steps, the 

functional relationships between the signals and the speed components 
are approximated by power-series expansions by the presumed-small 
components Uo/Ui, U 3 /U 1 and the small fluctuation quantities jui/Ui . 

As time averages we get in the first case a criterion for the posi- 
tion of the "flow-iine related" coordinate system and ±n the second 
case, with sufficiently numerous, sequentially set orientations of 
the probe, an equation system which can be solved for the desired 
average fluctuation quantities and the average speed U, 

if W 6 neglect enough higher terms* The condition to specify the 
average flow direction is thus recognized as an approximation, 
and it can apparently lead to a wrong calibration of the system in 
higher turbulence. Depending on the order to which the series 
expansion is carried, one obtains various approximations in the 
determination of U and u] , of which the first (linear) and the 



second are best-known. In suitable systematic notation, the series 
expansion can easily be taken to the fourth order in order to derive 
special, higher approximations. 

Instead of beginning from the average values e, , 
and of the X-wire, one can also use the quadratic signals 

but neglect higher terms. This approximation 

and the conventional approximations of first and second order, 
as well as a special fourth approximation of Vagt, are subjected to 
a numeric test since an experimental check of the various end- 
formulas is not possible. No better measurement method is available 
for this. As a great simplification, in a computer program, stochas- 
tic speed fluctuations of a turbulent flow are replaced by periodic 
fluctuations describable by step functions which are then easy to 
integrate, whereby the turbulence intensity varies and 

through suitable parameters, the quantities and 

u[u^/ (V assigned to different, fixed values. 

Even the convexity factor |T7][^/ap can assume different values. 

The test leads to the following results: All studie d an alyses 
lead to a systematic error in all turbulence quantities iujuyu* at 
high turbulence levels. This error is as a rule, greater in the 
analysis of quadratic signals than in the analysis of standard 
signals, so that the former can be discarded. The error is linked 
with the approximation character of the criterion for finding the /90 

direction of the average speed. The turbulence intensities 
are as a rule, computed with a greater accuracy than the shear 
stress terms lU , if these are different from zero. The 

different approximations are affected with errors of comparable 
size for large turbulence. In particular, it is not true that 
the higher approximation is generally the better one. For the 
limited range of turbulence intensity of 0-0.05 (turbulence degree 
0-227o), the max. error in the determination of turbulence intensi- 
ties tip/u^is about 30%, whereas the error in the shear stress terms, 
can be up to 100%. The deviations from the exact values depend on 
the turbulence structure without any discernable trend. Thus , no 
method can be given for correcting this error. However, from the 
calculations the maximum error can be found for each turbulence 
intensity. 

For minor deviations in the main direction of flow from the 
direction of the probe axis of an X-wire probe, the usual analysis 
of the hot-wire signals can be slightly modified and gives — with 
various-order approximations — the direction of the average flow 
and the turbulence quantities. For larger deviations, alignment 
of the probe is required so that when measuring a flow field, all 
measured quantities have to be transformed if they are to be related 
to a fixed coordinate system. These transformations can be given 
in a generalized form. 

Besides the average values 01“^ of the speed fluctuations, we 
are naturally also interested in the corresponding averages of 
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individual frequency components, or the spectral intensity densities 
17]'^ and the cross-spectra lU. As long as linear relations 

exist between the speed fluctuations and certain hot-wire signals, 
e.g. in the precise analysis via 3-wire probes or in the first 
approximation of conventional hot-wire anemometers with X-wire 
probes, the equations derived there apply equally for individual 
frequency components. Also, the equations of the second approxima- 
tion remain valid for frequency components and the unfiltered sig- /91 

nals are to be used only for the correction factors. 

Instead of X-wire probes, one-wire probes rotating slowly 
about the probe axis can be used if the main flow direction is 
known approximately. For the signal analysis however, low turbulence 
intensity must again be presumed, since a series expansion by the 
fluctuation quantities is unavoidable. The max. information, 
namely all quantities U.m = 1,2,3) is provided by a probe with 

slant-set hot-wire whose axis lies exactly or approximately per- 
pendicular to the main flow direction. With this method too, an 
error is expected in the turbulence quantities at higher turbulence 
levels due to the series expansion; this error is of the same mag- 
nitude as for the X-wire probes. 
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Tables 


Table 1: ERRORS ^,, AT os AND ISOTROPIC TURBULENCE 
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Table 3: ERROR AT AND 3-DIMENSIONAL BOUNDARY LAYER 



Figure 1: Normal One-wire Probe with Speed Components 
and Instantaneous Setting Angle 

Key: 1-speed components 2-hot-wire probe 3-probe shaft 4-prongs 

5-mantled ends of the hot-wire 6-effective part of the hot wire 


Fig. 2: X-Wire Probe with normal 
Prong Leads and Typical 
Dimensions . 

Key: 1-probe shaft 2-prongs 
3-dimensions in mm 
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Fig. 3: Orthogonal Three-Wire Probe ("Triangular Probe") 

Key: 1-front view 2-side view 3-top view 4-isometric spatial 
illustration (reduced) 



Fig. 4: Orthogonal Three-Wire Probe ( "Tripod" ) . 

Key: 1-front view 2-side view 3-top view 4-isometric spatial 
illustration (reduced) 
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Fig. 5: Plane Three-Wire Probe ("Right angle triangular probe") 

Application 1 

Key: 1-front view 2-side view 3-top view 4-isometric, spatial 
illustration (reduced) 



Fig. 6: Plane Three-Wire Probe ("Right angle triangular Probe") 
Application 2 

Key: 1-front view 2-side view 3-top view 4-isomet:ric , spatial 
illustration (reduced) 
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Fig. 
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Correction Factor F for Various Turbulence Intensities and 
Anisotropy Factors 
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Fig. 8: Correction Factor for Various Turbulence Intensities and 

Anisotropy Factors 



Fig. 9: The Corrected, Average Speed 
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/log 


Fig. 


Key : 


Fig. 


Key: 


10: Turbulence Quantities in an Airplane Wing Boundary Layer, 
Hot wire measurements by Elsenaar and Boelsma [33]; 

_y= local coordinate perpendicular to the surface 
= boundary layer density 


1-two-dimensional turbulent boundary layer 2-three-dimensional 
turbulent boundary layer 3— still predominantly two-dimensional 


turbulent boundary layer 




A C, • U^/Up 

o Cj t u^/u^ 
o C,j « u'j u'j/l 
A C,3 » Uj U*j/{ 

B Cjj • u '2 u*j/( iu^Vi/y^) 


11: Turbulence Quantities in a Boundary Layer in Front of a 

Cylinder Standing on a Flat Plate, Hot-wire measurements 
by Dechow [35]; (y = local coordiante perpendicular to 
the plate ) . 


1-two-dimensional turbulent boundary layer 2-three-dimensional 
turbulent boundary layer 3-still predominately two-dimensional 
boundary layer 
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Fig. 12: Turbulence Quantitiec in a Boundary Layer Generated by a 
Trip Wire on a Flat Plate, Hot wire Measurements by 
Charnay [34]; (y = local coord, perpendicular to the plate 
boundary layer thickness) 

Key: 1-two-dimensional turbulent boundary layer 



Fig. 13: Turbulence Quantities in a Shear Layer of a Fully Expanded 
Axial-symmetric Free Jet, LDA Measurements by Lehmann [36] 

Key: 1-axial-symmetric free jet 2-intermittence 
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Fig. 14: Relative Speed Fluctuations (i=i.2,3)for Isotropic 

Turbulence; analysis of "normal" signals: □ < < o.?s 



Fig. 15: Relative Speed Fluctuationsup’/u^ U = i,2.3)for Isotropic 

Turbulence; analysis of "squared" signals: o < < 0.25 


Key: 1-isotropic turbulence 2-"normal" signal 3-"squared" signal 
4-linear approximation 5-second approx. 6-fourth approx. 
7-approximation 




Fig. 16: Relative Speed Fluctuations (i = i , 2 , 3 ) for Isotropic 

Turbulence; analysis of "normal" signals; o < uj^/u’'< n .05 



Fig. 17: Relative Speed Fluctuations . 2 . 3 ) for Isotropic Tur 

bulence; analysis of "squared" signals s < o.os 

Key: 1-isotropic turbulence 2-"normal" signal 3-"squared" signal^ 
4-linear approx. 5-second approx. 6-fpui:th approx, /-approxi- 
mation 
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oos 


Fig. 18: 



Relative Speed Fluctuations (i=i, 2 ) 

Boundary Layer; analysis of "normal" 
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for two-dimensional 
signals; s < u .25 


Fig. 

Key: 



l_two-dimensional boundary layer 2-"normal" signal 
3-"squared" signal 4-linear approx. 5-second approx. 
6-fourth approx. 7-approximation 






Fig. 


Fig. 
Key : 


20 : 



Relative Speed Correlation v|ui/u* and Rel. 

2D boundary layer; analysis 
''squared” signals ;j osirp’/u^sD.ps 


/HO 


origimal page 'is 

OF POOR QUALITY 


Difference 
of "normal" 



1-two-diminsional boundary layer 2-three-dim. boundary layer 
3-"normal" signal 4-linear approx. 5-second approx. 6-fourth 
approx. 7-"squared" signal 8-approximation 
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Key: 1- also 4th approx. 2-no fourth approx. 3-Three-D boundary 
layer 4-"norTnal” signal 5-linear approx. 6-second approx, 
7-fourth approx. 8-approximation 
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Fig. 24: Relative Speed Fluctuations for Various Tangential 

Cooling Factors k and various Convexity Factors^r ; 
analysis of "normal*' signals: o s ITp'/u^ < o.os 


Key: 1-isotropic turbulence 
4-approximation 


2-*'normal" signal 


3-linear approximation 
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Appendix A: Linear Equation Systems, Matrix Multiplication, Deriva- 
tion of Equations (5.1) to (5.4) 

A linear system ot equations: 

^ a.K Xk - ^ ' 

is abbreviated as . , 


Here, A represents the matrix j x represents the matrix 

(x, )k=i,K, 5 b represents the matrix (b.)i=i,i. . 

Written explicitly, with K=3, 1=4 for example: 



3x1 

/ 


"I 

X, 


1 1 



3zj 

O 

Xi 


J 









The o multiplication is thus performed so that the individual 
rows of A are multiplied with column x like vectors handled as 
scalars. “ In this notation, linear equation systems are easy to 
handle since the rows can be multiplied by any scalars (A 0) and 
added . 

The equation system (5) is written for the orientations 

oi=: 7 i/ 4 ; 0 = 0, it; it/2, -tt/ 2; ir/^, -3lt/4i -7r/4,3it/4 aS folloWS: 


/114 
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0^4U,'J 

^ * »Z )( Uj + III ) 


I <(-; +w/;^ 

I C/j J. )^ 

(C/3+ up^ 

( Ur-Vnti+u,7 

(Uzi-uliiUs + up 


See (6.1) to 
(6.4) 



From this follows 


f 


2 C e| - e/) 


^ + €*) 

iCe,^~eP 

5 ( + e-,*" + e^V e j.® j 

5 4 - Cf - cj' — 

f C - e^ + e/_ 4^; 

Xie^ ~e^ ~ e/ + 

* j 


M 

M 

VJ 1 



liiA 




Appendix B: Proof that System of Equations (5) cannot be solved formally 
for the Six Unknowns (l,m=l ,25 3), when “= 


The question is whether there are six directions (<j,e j), 
so that: 


. c.fft • • j,. c,j-< .r,-K S,Hd. J,- , -I 


=+ 0 


Without affecting the universality, we can set e^-a . Then we have 


/•I,* 1 . * , !<•»»» ■ (. /t f ( Cot *o( ^ ■ 

2)fr = COJ.C .f.Vret • (fc 

1 ^ y O ^ 

A 4- Ck^-'T} €< ' 

(co/e, 6>v x/. €>.• , x,%. t J 

. ,'^v^ .<.k^-y,)t.4,^k ■ ck*-»;x/-v • 


</e^ 


A 

( o 


Ck*-^)ccJ« x<"« • (k*-c;cot« j.„« • 

Cot Cots’.- x.« , X-.V. S’/ U‘ o 3^.-v ^ 


• »Je4 ^ C ^ 


since the first and second columns of the last matrix are identical 
except for the sign, thus they are linearly dependent. 
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Appendix C: Calculation of Coefficients 
(7.1) to (7.6) for k=0. 


(«i*6i)of equations 


/117 


If we select the orientations (a.e) = (|, o), (l. „), (i, f), (§,-|), 

7^’ ^3’ ■ 7“^ and denote the corresponsing signals as e^ to e^, then 
for k=0 we have : 









•«l 


f \ 

4 

H 

i 

3 

V 

A 

4 

jf 

^ VI 
2 

il 

0 

0 

0 

0 


ui- 



V 


A 

2 







-V 

7 

f 

■? 

s 

•/J 

z-a. 

■A 

A/ 

■fj 

o 



2 

J. 

li 

s 

! 

s 

g 

■n 

2-H 

3 

> 

■fz 

2(Z 




Cy 

A 

V 

7 

S 

7 

? 

•/T 

~2-G 

/I 

- V 

•fi 


Uj,Uj 


e/ 

3 

s 


■n 

3 

{L 


UjM, 



H 

\ 

i 


2-fz 

V 

zfz 





By addition and subtraction of rows, we obtain: 


r > 

yf O 0 o o 


✓ N 

o 

o 

o 

<N 


u/ 

0 A —A 0 0 o 



0 0 0 4 0 0 

O 

■2ia^H2 

0 0 0 0 A 0 


2 

0 0 0 0 0 A 


2 W3L(« 

^ / 


V .. 


1 [, + 

[ - 2 e-1^ - ^ ®3 + -t J 

oT-^TVI f - 6 - e/ ^ 3 ] 

[ et + - e-s - ei ] 

e| - fiv - -3 e-/ e 4 ] 


/118 


76 



and furthermore 
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Appendix D: Hot Wire Analysis without Root Expansion, 

and Handling of Equation (42) 

Let U2=Uo=0 and presume u>>|u|| . From eq. (27) it 
without any '^approximation: 


4. 



// 


i ^ 

/ 

3^a ^ 2 

QAA ^4 






ifi' 4 

^ i. £5/14 

>2 



f wi® + 

Ml? ^ 


Hi- + 



r ^ 

3><2. 

2 

-f Z 

iaiii ■*■ —nr 

SZ. U2 " 


and in consistent expansion out to the second order: 






, . ui* ^ ±1^ i 

o ^ 4- 2 'TT^ ^ 

/I ^ 2 -jfz ^ u 


/Z 


^^4 \J4 


B-4/Z. 









3^z 

u>al 

4- 

/(2 

£V<'t 





!S 

+ 

xl2 


(/<'■ 


P ^ ^ ^ «5** 


t»/ ) I 


^ 4 6 -f C2 


a» _ ^ (3^2 f ) “iYa +• 


\' 


From this follows equation (42): 


Derivation /H9 


follows 
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^ i -2 3-fw ^ *2 <3<^ ^3^ ^ '“4' a.f2 ^ 

^ ^a>ij ^ (2 34i,a2a •^'4'a^J)^ ^2 a,(.fa33+ 4'a4), "'^■^ a^.f a^j. ^ 
63/1^ (2c?^32i — 4'£?^5) I 2 3^ 3ii — ^ Q^t) 0 

f / * ' / 



and furthermore : 
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/121 
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and thus : 



7122 
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There is a very great formal similarity of this equation 
system with equation system (32). , For the usual eight orientations 
of the hot-wires, the following equations result: 






u'^V} 


(Jj- 






A 

lA+k^i^ 


r {.cUc^)^ 


'I T 

i ^ 

-S' J 

(A-t-k‘)^ 1 

lA4k‘ 


fc, - ?i)' 

t yf-hk^ \2. 
VA- kV 

1 * 

Cc^-Z^)' 

•f 

(«Sj- 


/ 


-I 


d. ( 1-^. {cJ_-Cc,‘f- 4- CCf-cJ^- fc^'- Ct)^ -C<r>-c>;^ 

+ A-kW ~~f — — ^ 


(A+k^) 

_d / \ - c,'*) 


L^T( A- 2 (4^ U/ H- 2 =- 

“ tC e-3'..gr,^; ^ 

2 ^ o/j + 2 + 

^ / — 2 . 2 

a C *^5/^ + 


- ^ yj-t-k^ dzAi A t -r . — 2-1 

+ ^ a -V- • u, " a C j 
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Appendix E: Properties of the Step Function g(t',s,e) 


g(t ' , s,c) = J 


" V / 0 ^ i' & s ^ 

, s ^ i' er 


^/<r 


/d _^r 

y s - e 


^ G~ ^ i ^ 


continues periodically on the t'-axis 


O S ^ (S' ^ 4 


3 • J oU' 


3) 


s- 


/"w T’ 

r/rry ^s{fr^' 




s 


— ^ tt _ SL 
to J C5" 


7T 


J S' 


i) 3* ° ^ J+ fe=^C(T- 




» J ( -f - d:) 4 — . «. , £f T- ^ 

-* or S - ^ 

— + w yi 

a" ~ ^f-'k) • <■«>■-•"■> 




~ 6 " / ( — 3 /e- . 


4) Assuming that 2s £ 6 1/2 


a) for D < r < s: 


**jr yt T ( + at> ^ 
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b) for s < r < <S-s: 




— s -h (G'-s-r) i>^ “ 


c ) for 


^•s< r •< ^ 


g(i', J-, «■> • g( t'-f, - — ( S'- H • dfe - 

iir ^ — f 4 C6‘-r>at + C«r-J-) fc^J . 


d ) for <r < 


i •• 


dli',S,G) ■ 


e) 


3 ir~r^^,G) =r JU^T7rjJT~JTp;7JJ 

“ 3^',s,<5) ■ 3ir+r,s, 6-) . 


Thus, together we have: 
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/f ~- I Irl ( + e>b * b^)], 0 ^/f/^ J' ; 

/f - [ S ( ^L) -Hrl^h^] ^ J 4, tr! e'-i; 

^ _ [ jr - (<S--trl) Ob 4 - 

+ ( <f - i-) b D , £■•■ 4 ^ Itrl i- s'- 

0 


and continues periodically 
If we substitute in a,b, then we obtain: 




A - 

/ r J. \ 0 ^ tr! *■ s ; 

4/X - '/ 

Irt 

V ^ / r/ ^ 


/.r ^ 

- 4 

1 r/. s ^ If^f ^ ^ ^ 

<o ' 

O 

S' ^ t ; 

/ 


and continues periodically. 
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Appendix F: Computer Program to Check the Conventional Analysis 
and the Analysis without Root Expansion 

Nomenclature: 

Name in com- normal quantity 

puter program designation 


K 

k 

S 

s 

T 

e 

Gftl’l 

r 

C2 

Cj 

C3 

Cs 

SA1 

sign( Ai ) 

Sfl2 

sign( Aj ) 

SA3 

sign(A} ) 

C(1 ) 

C j 2 

C(2) 

Cjj 

C(3) 

C 2 i 

U ( 1 , N ) 

(ei -ij )/tei +Ej ) 

U(2,N) 

(ia-ii, )/(ej+i4 ) 

U{3,N) 

(es-Gfe ) 

Q(A,N) 

(B7-e> )/<e7 + i6 ) 

LAl<)( ) 


MY ( ) 

V 

NY( ) 

V 

W1 

irp-/u^ 

W2 


U3 


UA 


U5 

ITJTrj/u^ 

W6 


W7 

t?JUT-/u2 

) 


E2 

^1/2* 

EPS1 





EPS2 

(c'-c^')^ + d77^ 



F2Q( ) 

F * 

0 
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MNF. 


1. 

6000003 

c 

c 

c 

lAnalysis of hot wire signals of turbulent flows 
Type and dimension Instructions 

INTEbER 0 


coanooB 


REAL K»LA3< 11> > 47 t 11» , HY< iH , K4P,K2 , LA40 (11! »«7Q ( 111 , NrQlllI 

3. 

OOOOOCB 


Ji«ENSl3N R( 3 ! , AL ( 8 1 , C ( 3! , AU K < a • 3 ) , G ( 3 . 17 ! , 5M ! 8! , 



c 

c 

1 GS(4!5&D*4),Q(4,ll!,4A(5,17J,EP(7>,Wl7,a,ll»jF2!ll)» 

2 f E'( 8} ,F < 3 ? 171 ( X ( 16! . 9x { 17 1 > X4< 17! , GQ( 8,17) .F2Q tin 

(Input of parameters and calculations of signals 
1 Isotropic turbulence (beginning ) 

*» • 

0021313 


r=i./'. 

5. 

00 4 64 


R!ll =T 

6. 

0046408 


S. < 2 i = C . 

r. 

B14647B 


Rt3»=2.»T 

8. 

804650B 


UO 25 LS=1,2 

9. 

Cui4 66 ^ B 


i=l./b.-FuOAT(LS-ll*tl./o. -4. 1.6531 

iC. 

0046569 


4lk12=1,C 

It. 

004 65 73 


AlP13=1.C 

12. 

004 66 

c 

i Isotropic turbulence lena; i. 

13. 

0046613 


SA2=SIGMtl.',Ai.P121’SAl 

1»». 

0046636 


SA ? = $IG ►((!., 4lP13)*?A1 

15. 

0 0 4 6668 


Oq3L— 1^3 

lb. 

0046719 


X(i_) =6-Rt L) 

17. 

Cu46723 


X { J.i. ) = T-r( L) 

la. 

U046758 


X (a+L ) =1. -6 (U< 

19. 

0 J 4 7 0 r S 



2b. 

u 0 4 ? B 


3 < tl2+U =1 . fT-R (U 

21. 

CO 4 7 0 79 


X ( L61 =r ,c 

22. 

D0470 73 


X 1 17) =1.0 

23. 

504^100 


X ( 18) =1.0 

24. 

ncu7iio 


ArtONlAfldi 

25. 

t* 3 4 7 m t3 


JO lb -'=1,18 

26. 

bC<.715b 


i F (X t N) i . 1 . 0 ) X ( 4) = , C 

27. 

fC4721B 


aF (X ( >'l .61 . 1 . j ! X ( v) = 1 

23. 

C:h7243 


to 3 i N'JE 

29. 

CC47263 


J4 15 M=l,17 

30 . 

0447303 


9 X(«) -X (4 H> -X ( 41 

31. 

0447113 


15 XX' (H! = {X(4+1! 7X (M( ) /7 , 

32. 

00 47353 


DO 25 L<=1,2 

33. 

bO>.77oB 


K = u .E^ffAT «lK-U 

34. 

CJ474C3 


0a = Ai.P12»»2 

35. 

£3^.7418 


Ci=ALPi3»x2 

36. 

0047^^>» 


=1=4. *ATAN{1 .) 

37. 

(.34714^-3 


DO 1 1=1, 8 

3«. 

04‘.7513 


1 Ai. tl) =Pi/4. 

39. 

C047548 


Tt !1! =0,0 

40. 

004754B 


rF(2!=Pi 

41. 

0047553 


TE (3 ) =P;72. 

42. 

0047573 


Ttu t4) =-PI/2. 

47. 

00476I:B 


TE(5)=-:/4. 

44. 

CJ4761R 


TE (.6) =-o. »PI74 . 

45. 

Ou 4 7 63B 


r£(7)=-p:/4. 

6* 

C0475<.b 


Tc t3)=3.*Pi74, 

47. 

C047668 


GAM= t / (S» (i-i) ) -3 ./T 

43, 

0047713 


00 4 1=1,6 

49. 

0047733 


4=aL(l) 

5u « 

004773B 


TH=TE (X! 

51, 

004775B 


■<2 = K’‘*,:,-1 , 

52. 

9 3 3 (J Q C 8 


AlHU ,IS = 1.4K2*(C0i (A? )»*2 

53. 

0050073 


Ai..4(c,I)=l.'X<2*tiIN(A)«CGS(THn»*2 

54. 

005017B 


AL4(3,I)=i.'>'K2'»{jI.4(a)‘S:suIhM'x*'2 

55. 

005027B 


AL4(t,i) = !<2»CJi<A)»bXN(A)»C0StTH) 

56. 

r 0 5 0 4 0 3 


AlM( 5 ,I! = '<2»CL/i(«) 'il,, (A) »sin(T-i| 

57, 

005 05 23 


Au-!{fc,i! = <2»(SX'4(A) ) *»2'C0S(TM)*S:'4(TH( 

58. 

0050548 


4 0 u.»r I NUE 
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59. 

0050703 


jO 2 4=1,5 

60. 

1050728 

2 

ca.) = c.u 

61. 

0050743 


UU 5 4=1,17 

62. 

0050768 


Oa 55 L = l.,5 

63. 

0051403 


FiL,.Ml = FF(S,T,PlLi,X’'(''n 

64 . 

0051003 

55 

CCi7lNU£ 

65. 

C05112B 


C{1')=C(1J+DX{M)*F(1,m)>!:{2,M) 

66. 

0051178 


G<2i=C<2l +0X(H»*r d,M) *F<3,MI 

67. 

005123B 

5 

C<3)=C(3)»DX<M)»F(Z,*')*Ft3,«) 

6S. 

(>451333 


0 d) =SI8N (1 ALP12J »C( 1) 

69. 

0051768 


u(2)=SIGN(l.,AL'>131'C<2> 

70. 

0J514CB 


3(5) =8lGl4<l . ,4LP12 1 »Bi.GN( 1 . ,Ac»l 31 ■»(; t 3) 

71, 

OO51450 


OU 2C N=l,il 

72. 

0051478 


41 = ';Ai*5Q<T 10, ,25*FuUAT (H-l> 1 

73. 

0051558 


«2=AuPlt‘Al 

74. 

0051578 


Ao=AlP13*A1 

75. 

0o516C8 


80 6 1=1,8 

76. 

0u51623 


JG 6 M=l, 17 

77. 

0051658 


G4=A..Md,I)*<l.*Al«=d,4)»*^7-ALM(2,:>* 
( ■»,!) * ( A3*F (3 .9) ) "*2 
2 + 2.»HL'r(4,:i*»l.»Ai’Pd,.!j»»A2*''(2,'() 
5<-2 .*Alm( 5 , I » » ( 1. + Ai*P d ,•:)) ‘Ap*' ( 3 , .-1) 
4f2,*AL-((6,i.)*A2"'<2,M) *A5*Fd, •'1 

78. 

''052148 


■5Q(I,M)=;-2 

79. 

(.0521 73 

G 

6 

Gd , •M = S0PT » 

Calculation of time history 

.80. 

345232b 


0 U 7 1 = 1 , 5 

81. 

0052o33 


^ = : • O' 

82. 

CO 5 2 p38 


8 m=1,17 

83. 

0052568 

Q 

“=•'^ uxr<) *G ( : , '<) 

84. 

u052458 

7 

G.-.( I 1 =^ 

65. 

0o5?52o 


° J-1,4 

86. 

''•J575 03 


GO ( J) = (^M (4* J-1 ) -v»4 <2* J) ) / 2. 

87. 

C052569 


( Jl = (uM (B^ 1) +G''U2’‘ J) ) 7 2. 

8d. 

10520^8 

i 

; ( J, I =G'J( J) 7G J ( Ji 

39. 

C-15264E 


L. A . ( N ) = G6 1 2 ) 7 i P (1 ( 

90. 

10 5 2r> 7 8 


tx ( f! i -=GS ( 5) / jS ( 1) 

91. 

4 ^ 5 J •* ’. 8 


.X (N> {..)/■. 5 d> 

92. 

205 27 58 


w’G 1'4 1=1,8 

9 5. 

0d52753 


Uu lu M=l,l7 

94. 

0352779 

1 

2 i,A(X,Mi=4(A,*'0— ^'■1(.) 

95. 

C 0 5 3 0 7B 


jn 11 j=i,7 


96. 

gr, 

98. 

99. 
108 . 
101 . 
102 . 

103. 

104. 

105. 

106. 
107. 
188 . 

109. 

110 . 
111 . 
112 . 

113. 

114. 

115. 

116. 
117. 
IIS. 
119. 
120 . 
121 . 
122 . 
123. 


0J5311B 
0j57i33 
QJ53158 
005322b 
025 3246 
C05331d 
C35335B 
005.5419 

0053458 
045 2528 
QC5355B 
OU5''6rti 
0 053648 
0053708 
0053748 
0054UCB 
00540 48 
0u54O«B 
0054118 
00541 68 

0054238 

0Q5427d 

0054328 

0054358 

0054438 

0054503 

0054569 

0054618 


€Ptj>=).: 

00 12 *^=1.17 

iP (l»=t8(U*t0A<2)H)*-0Mljd)l**2 * 

cf (2)=rP(2)»(u7t2j!-H-0A<l,''>)’»2 * 

,:4(3}=irp(T) + ((;A(4.w|-GAt3.Hn**2 * 

.■0(Z!=-f-(4)fl54(2,'-')**2-«fl<l.'it*»2) 
£P ( 5 ) =cP (5) KGM(4,Mj»»2-iA<3,M)»*2) 

£Pt6J=56t6?Kl>M(6»Ml-GA(y,MJ»*»2 * 


12 


El- (7 1 =pp ( 7 ) * (G A ( d 1 :•!) - GA ( 7 

, X 1 ) • ' 2 

fcl=tPdi/(d0u.*G'd)”‘2) 

• 2 : 0 . 

c2 = £P (2)/(8d:.*GSd)’‘*2) 

»20 2. 

=.5 = £P(3!7(d0 0.*G5(2)»»2> 

•2 3i.- . 

£A = F.P (4) 7 ( 8 00. '■GS( it »♦ i) 

•200. 

£5=FP!5!7(300.*GS(2)**<) 

'►20C . 
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Appendix G: The Second Approximation for the Hot Wire not Aligned /I36 
in the Main Direction of Flow 


According to ( 27 ), as a second approximation, we have: 
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With the usual X-wire orientations, we thus obtain a linear 
equation system of 20 equations for the (formal) 12 unknowns 
2 ,U,U2.U2U3 .UjU^ , tjy^, u^, i[j|uJ, uJu^, . 

From these equations follows: 
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and other equations which will not be needed. From the given 
equations, the final result can be derived: 


C ^ 2. 
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The formulas for the 1st approximation follow formally from the 
2nd approximation, if we set: 

Fj = 1 and = Ss-/c ■= s^/e 

d 


It is easy to see that for the case when U 2 =Uo or ‘^i/2='^3/4®‘^s/6‘''^7/b-''° 
the formulas (G2.1) to (G2.6) and (G2.9) are transformed into 
the old formulas (33.1) to (33.6) and (36) of section (4.1). 

On the other hand, we have: 
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for the transition from turbulent to laminar flow. 
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Unfortunately, the formulas are very cumbersome. In many cases, 
they can be made much simpler; If we can assume U^=0 for symmetry 
reasons and by neglecting tangential cooling (k=07, then we have 
for example: 
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We see that changes compared to the aligned X-wire occur only 
in the equations for iJlp*’, ZJ^, ZJ!]J and F^, where the correction 

terms for IT!p:, Zp' are of fourth order, and are thus negligfble. 
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